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SUMMARY 


Analytic  approximations  for  many  aspects  of  airblast  from  nuclear 
explosions  are  presented  in  this  collection.  Employing  shock  over¬ 
pressure  as  an  independent  variable,  shock  front  parameters  are  fit, 
including  shock  density,  the  specific  heat  ratio  for  shocks  in  sea- 
level  air,  shock  velocity,  peak  particle  velocity,  peak  dynamic  pres¬ 
sure,  shock  temperature,  and  normal  shock  reflection  factors  in  sea- 
level  air. 

Free-air  (and  surface)  burst  nuclear  blast  parameters  are  approx¬ 
imated,  including  such  relations  as  peak  overpressure  versus  range; 
shock  radius  versus  peak  overpressure ;  peak  overpressure  versus  time 
of  arrival;  time  of  arrival  versus  peak  overpressure ;  time  of  arrival 
versus  shock  radius;  positive  overpressure  duration  versus  time  of 
arrival;  peak  overpressure  or  range;  overpressure  impulse  versus  peak 
overpressure;  overpressure  and  dynamic  pressure  versus  time  and  range; 
overpressure  impulse  versus  peak  overpressure;  overpressure  and 
dynamic  pressure  versus  time  and  range;  overpressure  impulse  and 
dynamic  pressure  impulse  versus  overpressure  or  range;  maximum 
(fireball)  temperature  versus  peak  overpressure ;  and  negative  under¬ 
pressure  versus  peak  overpressure  and  time. 

Some  blast  parameters  as  a  function  of  burst  height  are  approxi¬ 
mated,  including  overpressure  and  dynamic  pressure  versus  burst 
height,  ground  range,  and  time.  These  fits  are  compared  with  results 
of  detailed  computer  calculations. 

The  equations  of  state  for  air  are  fit,  including  both  the  ther¬ 
mal  (P  =  RTp )  and  caloric  [Ep  =  P/(Y  -  1)]  relations.  Some  of  the 
fireball  thermal  radiation  characteristics  are  given  in  further 
analytic  approximations  for  surface  bursts,  air  bursts,  and  high- 
altitude  bursts. 
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PREFACE 


This  collection  of  airblast  fits  was  prepared  at  the  urging  of 
Dr.  George  W.  Ullrich  (Defense  Nuclear  Agency/Shock  Physics  Aerospace 
Systems),  and  Mr.  Clifton  B.  McFarland  (Defense  Nuclear  Agency/Shock 
Physics  Strategic  Structures),  who  perceived  a  need  for  a  single 
document  summarizing  the  available  analytical  approximations. 

In  compiling  these  fits,  the  author  relied  upon  many  sources  and 
individuals.  In  particular,  the  years  of  detailed  curve  fitting  so 
diligently  pursued  by  Stephen  J.  Speicher  (a  coauthor  on  numerous 
previous  reports)  are  liberally  represented  in  the  approximations 
offered  in  this  document,  and  in  the  reports  noted  in  the  reference 
sect  ion. 

Special  thanks  are  due  Gilbert  C.  Binninger  of  Science  Applica¬ 
tions,  Inc.,  La  Jolla,  California,  and  Fred  Sauer  of  California  Re¬ 
search  &  Technology,  Inc.,  Pleasanton,  California,  for  their  thorough 
review  and  numerous  suggestions. 
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CONVERSION  TABLE 


Conversion  factors  for  U.S.  customary  to  metric  units  of  measurement. 


To  Convert  From 

To 

Multiply  By 

Pounds  per  square  inch  (psi) 

Kilopascals 

( kPa) 

6.8951 

Bars 

Kilopascals 

( kPa) 

100.00 

Calories  (cal) 

Joules  (J) 

4.183 

Feet  (ft) 

Meters  (m) 

0.3048 

Inches  (in.) 

Meters  (m) 

0.0254 

Kilotons  (KT) 

Terajoules 

(TJ) 

4.183 

Megatons  (MT) 

Teraj  oules 

(TJ) 

4183 

Kilofeet  per  second  (kft/s) 

Meters  per 

second  (m/s) 

304.8 

Kilofeet  (kft) 

Meters  (m) 

304.8 

Kilofeet  per  cube-root 

Meters  per 

terajoule 

189.2 

kiloton  (kft/KT1 /3) 

(m/TJ) 
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SECTION  1 
INTRODUCTION 


Even  before  the  first  nuclear  burst,  the  effect  of  height  of 
burst  (HOB)  on  the  blast  wave  was  considered  important.  The  Hiroshima 
and  Nagasaki  bombs  were  detonated  at  considerable  altitude  in  order  to 
increase  the  range  at  which  significant  blast  damage  might  occur. 

Blast  measurements  were  made  on  more  than  90  atmospheric  nuclear 
tests  in  the  period  1951-1963.*  A  number  of  detailed  hydrodynamic 
and  radiat ion-hydrodynamic  calculations  have  helped  to  explore  the 
physical  processes  involved.  Further  measurements  and  understanding 
have  come  from  simulated  or  analogous  blasts  created  with  chemical 
explosives.  From  each  of  those  sources,  we  have  gained  information 
that  has  required  further  interpretation  and  simplification  for  appli¬ 
cation  to  targeting  or  to  the  design  of  survivable  structures  and 
military  systems. 

As  a  consequence,  frequent  efforts  have  been  made  to  find  useful 
analogs  and  algebraic  formulae  to  represent  the  principal  blast 
features.  Some  are  derived  from  theoretical  principles,  some  are  the 
result  of  empirical  fits  to  observations  and  test  data.  All  attempt 
to  simplify  and  quantify  the  complex  physics  of  blast  waves. 

The  analytic  expressions  for  blast  waves  presented  in  this  report 
are  divided  into  three  general  classes.  The  first  class  (in  Sec.  2) 
includes  relations  that  can  be  expressed  for  shock-front  variables, 
independent  of  yield  and  burst  geometry.  These  Rankine-Hugoniot ,  or 
shock-front  relations,  become  more  complex  for  strong  shocks  in  air, 
since  the  air  molecules  begin  to  disassociate  and  the  air  atoms  begin 
to  ionize  as  the  strength  of  the  shock  grows  and  as  the  air  tempera¬ 
ture  and  pressure  rise.  Section  2  also  includes  the  appropriate 
equation  of  state  for  air  at  the  shock  for  both  caloric  and  thermal 


*A  brief  summary  of  atmospheric  tests  and  a  list  of  references 
are  presented  in  Brode  [1979].  An  analysis  of  the  resulting  peak 
overpressure  HOB  data  is  provided  in  Brode  [1981]. 
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properties,  i.e.,  for  both  internal  energies  and  temperatures. 
Specifically,  Sec.  2  deals  with  shock  density,  shock  velocity,  par¬ 
ticle  velocity  at  the  shock  (shock  wind),  peak  dynamic  pressure,  shock 
temperature,  and  shock  reflection  factors.  A  tabular  summary  of  those 
variables  is  included,  along  with  peak  overpressure,  the  ratio  between 
pressure  and  energy  density,  i.e.,  the  effective  gamma-minus-one,  and 
the  gas  "constant"  for  the  thermal  equation  of  state,  i.e.,  the  ratio 
of  pressure  to  temperature  times  density. 

The  properties  of  a  nuclear  explosion  in  air,  unaffected  by  any 
surface  reflections  (a  "free-air"  burst),  have  been  described  in 
detail  by  numerical  radiation-hydrodynamic  calculations,  and  have  been 
used  to  describe  surface  bursts  (assuming  a  surface  burst  is  equiv¬ 
alent  to  a  free-air  burst  of  twice  the  yield).  Some  of  those  nuclear 
blast  wave  properties  are  described  by  simple  algebraic  fits  in 
Sec.  3.  The  relations  provided  include  peak  overpressure  versus  shock 
radius;  shock  radius  versus  peak  overpressure;  peak  overpressure 
versus  time  of  arrival  (TOA) ;  TOA  versus  peak  overpressure  or  shock 
radius;  duration  of  positive  overpressure  versus  TOA,  peak  overpres¬ 
sure,  or  shock  radius;  overpressure  impulse  (in  the  positive  phase) 
versus  peak  overpressure ;  overpressure  versus  time;  duration  of  posi¬ 
tive  dynamic  pressure  versus  peak  overpressure  or  shock  radius; 
dynamic  pressure  versus  time;  maximum  temperature  versus  peak  ovei — 
pressure;  and  time  to  maximum  temperature  versus  peak  overpr essur e . 

Tabular  listings  are  included  here  for  peak  overpressure ,  range, 
TOA,  positive  phase  duration  and  impulse  for  overpressure,  and  dynamic 
pressure  and  maximum  temperature  in  the  shock  or  fireball  (all  for  a 
1 -KT  free-air  burst). 

Section  A  provides  fits  for  peak  overpressure  versus  burst  height 
and  ground  range;  overpressure  versus  time,  burst  height,  and  ground 
range,  peak  dynamic  pressure  and  total  positive  dynamic  impulse  versus 
burst  height  and  ground  range;  dynamic  pressure  versus  time,  burst 
height,  and  ground  range;  and  some  limits  for  dynamic  pressure 
decreases  (10  percent  less  than  surface  values)  as  functions  of  target 
height  (in  terms  of  ground  range  and  burst  height). 
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Section  5  provides  a  detailed  fit  to  the  caloric  equation  of 
state  for  air.  The  thermal  equation  of  state  for  air  is  included  in 
Sec .  2 ,  Eq .  ( 25 ) . 

The  principal  features  of  the  thermal  radiation  from  the  fireball 
are  given  quantitative  approximations  in  Sec.  6.  That  section  in¬ 
cludes  times,  peak  power,  energy  radiated,  and  atmospheric  transmis¬ 
sion  for  both  air  and  surface  bursts. 

Forty  years  of  blast  research  have  yielded  an  extensive  data 
base,  plus  considerable  confidence  in  many  of  the  simplified  analytic 
approximations  presented  here.  However,  airblast  research  in  recent 
years  has  been  proceeding  at  a  rapid  rate.  Consequently,  these  fits 
for  idealized  conditions  will  need  corrections,  improvements ,  and 
additions  within  a  few  years.  One  of  the  areas  being  investigated 
most  intensively  is  airblast  over  "nonideal"  or  "real”  responding 
surfaces.  There  is  already  a  need  for  simple  representations  of  blast 
waves  over  such  real  surfaces  for  use  in  targeting  and  survivability 
analyses.  The  simplifications  and  analytic  fits  offered  in  this 
report  may  provide  a  logical  point  of  departure  for  extensions  into 
nonideal  surface  complexities  and  structure  interaction  models. 

Appendix  A  provides  overpressure  values  as  a  function  of  time, 
range,  and  burst  height.  These  values  are  provided  as  test  cases  for 
those  who  have  programmed  the  overpressure  fit  and  wish  to  check  for 
errors. 

Appendix  B  is  a  Fortran  listing  of  a  program  to  calculate  the 
overpressure  as  a  function  of  time,  burst  height,  and  ground  range. 

While  accuracy  estimates  are  provided  wherever  possible,  there  is 
no  adequate  treatment  of  the  variability  and  degree  of  unpredict¬ 
ability  in  airblast  features.  Measurements  made  during  the  atmos¬ 
pheric  test  program  generally  supported  a  variation  exceeding  plus  or 
minus  20  percent  on  many  blast  features.  A  review  of  peak  overpres¬ 
sure  measurements  by  the  author  [Brode,  1981]  emphasizes  the  data 
scatter  and  the  lack  of  reproducibili ty  in  the  nuclear  blast  over  real 
surfaces.  Clearly,  the  analytic  expressions  presented  here  cannot 
reflect  such  large  uncertainties. 
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SECTION  2 

SHOCK-FRONT  VARIABLES 


For  shock  waves  in  an  ideal  gas  or  fluid,  all  of  the  dynamic  and 
thermodynamic  variables  at  the  shock  front  (density,  particle 
velocity,  shock  velocity,  dynamic  pressure,  temperature,  etc.)  can  be 
expressed  in  terms  of  preshock  conditions  (ambient  pressure  P0 ,  am¬ 
bient  air  density  p0,  ambient  air  specific  heat  ratio  Y0,  ambient  air 
temperature  0O)  and  a  single  variable  such  as  shock  strength  (e.g., 
peak  overpressure).  A  number  of  such  relations  are  derived  in  this 
section. 

At  high  temperatures  and  pressures,  air  is  not  an  ideal  gas, 
which  limits  the  utility  of  an  ideal  gas  formulation  for  strong  shocks 
in  air.  Since  the  effective  specific  heat  ratio  Y  changes  slowly  with 
shock  strength,  analytical  expressions  for  all  of  those  parameters  can 
be  found  which  follow  fairly  simple  forms;  some  are  not  unlike  their 
ideal-gas  counterparts.  In  addition,  the  influence  of  equation-of- 
state  changes  can  be  represented  by  a  fit  to  Y  at  the  shock  condition 
in  air. 

Fits  for  peak  dynamic  pressure  Qs ,  shock  velocity  Us,  peak  par¬ 
ticle  velocity  us,  peak  density  ps,  shock  temperature  6S,  and  normal 
reflection  factor  RF  are  offered  as  functions  of  the  peak  pressure  P3 
or  peak  overpressure  (APS  =  Ps  -  P0),  the  preshock  conditions  (P0,  p0, 
Y0,  and  ambient  air  sound  speed  C0) ,  and  the  shocked  air  specific  heat 
ratio  Ys. 

CONSERVATION  RELATIONS  AT  SHOCK  FRONT. 

Conservation  of  mass,  momentum,  and  specific  internal  energy 
across  a  shock  front  are  represented  in  Fig.  1  and  the  following 
equations. 
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Figure  1.  Shock-front  parameters  (in  shock-front 
frame  of  reference). 
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Equation  (3)  may  also  be  written  as 
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where  Es  =  shocked  specific  internal  energy, 
E0  =  ambient  specific  internal  energy. 


(A) 
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Implicit  in  Eqs.  (3)  and  (4)  is  the  ideal  gas  relation  for 
specific  internal  energy  E: 

E  =  P/[p(Y  -  1 )]  , 


and  (or) 

Y 

P/p  =  constant  ,  (5) 

for  adiabatic  flows,  in  which  Y  is  the  specific  heat  ratio, 

Y  =  Cp/Cv,  where  Cp  is  the  specific  heat  at  constant  pressure,  and  Cv 
is  the  specific  heat  at  constant  volume,  and  where  P  =  pressure  and 
p  =  density. 

SHOCK  DENSITY. 

Solving  Eqs.  (1),  (2),  and  (3)  for  the  shock  density  as  a  func¬ 
tion  of  the  peak  pressure  P3  leads  to 


Expressed  in  terms  of  shock  overpressure  (APS  =  Ps  -  p  ), 
the  relation  becomes 


(7) 


6 


This  relation  for  sea-level  air  shocks  is  illustrated  in  Fig.  2. 

A  fit  to  the  shock  density  ratio  ps/ p0  for  sea-level  air  shocks  versus 
shock  overpressure  is  also  illustrated  in  Fig.  2.  The  fit  is  as 
follows : 


ns  P 


86.06  +  51  39  tt  +  3  -48  Ott  ‘ 


0. 002985tt 


2.473 


—  —  1  T  '  ■  '  '  ■  — 

o  86.06  +  926. 3tt  +  243. 8tt2  +  2.325tt3  1  +  3 .57*1  0~8tt3* 221 


3  7  4  c  ‘ 


3966c 


1  +  71 .82;  1  +  1  ,136,830c 

0.004599c3,893 

1  +  2.055*10~5c6,1 39 
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9.486c 


2.567 


1  +  5.911c 


5.127 


(8) 


where  it  =  APS/1000, 

C  =  APS/106, 

and  APS  is  in  pounds  per  square  inch. 


The  fit  is  within  5  percent  for  all  values  of  the  overpressure. 
At  very  high  overpressures  (>200,000  psi)  the  adiabatic  shock  assump¬ 
tion  is  unlikely  to  be  valid.  For  strong  nuclear  blast  waves,  radia¬ 
tive  transport  in  the  fireball  can  dissipate  or  redistribute  shock 
energy  and  pressure,  thus  allowing  further  compression  and  higher 
densities. 


SPECIFIC  HEAT  RATIO  FOR  SHOCKS  IN  SEA-LEVEL  AIR. 

The  effective  Y  for  air  shocks  (in  standard  sea-level  air)  is  a 
slowly  varying  function  of  the  shock  strength.  That  dependence  is 
illustrated  in  Fig.  3-  For  sea-level  air,  Ys  ranges  between  1.17  and 
1.67.  It  starts  at  Y  =  1.4  for  weak  shocks,  drops  to  1.17  around 
10,000  psi  (11 ,00°C) ,  and  eventually  rises  to  1.67  at  very  high  pres¬ 
sures  0107  psi).  According  to  Eqs.  (6)  and  (7),  this  means  that  the 
shock  density  rises  to  nearly  13  times  ambient  at  around  1  kb 
(15,000  psi),  but  is  limited  to  only  4  times  ambient  at  extremely  high 
pressure  ( > 1 07  psi). 
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Figure  2.  Shock  density  versus  shock  overpressure  in  standard, 
sea-level  air  (p0  “  1.293  kg/m^). 


Gamma 


Figure  3.  Specific  heat  ratio  (Y  =  Cp/Cv)  for  shocks  in  standard, 
sea-level  air  versus  shock  pressure. 
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An  approximation  to  the  specific  heat  ratio  for  shocks  in  sea 
level  air  as  a  function  of  shock  strength  is  as  follows: 

Y=1+A+B+C+D+E+F+G,  (9) 

s 


where  A  = 


0.4  +  2 . 39 9 tt  +  1 1  6.4tt 
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3.742  ’ 


G  = 


-9  8 

3.931x10  c 

1  +  2.587x10_8c8 


and  where  it  =  APS/1000,  in  kilopounds  per  square  inch, 

C  =  tt/1000  =  APgxlO-^,  in  megapounds  per  square 
inch. 


This  expression  fits  to  within  3  percent  of  the  curve  in  Fig.  3, 

which,  in  turn,  is  based  on  a  fit  to  the  equation  of  state  for  air - 

good  almost  everywhere  to  less  than  5  percent  error  [Brode  and  Parkin, 
1963].  (See  also  Sec.  5  of  this  report.) 

As  mentioned,  however,  at  very  high  pressures,  temperatures  are 
also  high,  and  radiative  processes  allow  even  higher  shock 
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compressions.  In  a  flow  field  dominated  by  radiation  transport, 
shocks  can  reach  much  higher  compressions  because  radiation  can  carry 
away  the  shock  heating,  thus  lowering  the  shock  pressure  that  would 
otherwise  limit  compression.  For  that  reason,  these  adiabatic  shock 
formulations  may  be  inappropriate  for  nuclear  blast  waves  much  above 
200,000  psi  in  sea-level  air. 


SHOCK  VELOCITY. 

Similarly,  from  Eqs.  (1),  (2),  and  (3), 


U 

s 


C 

o 


Y  +  1 )AP 
s _ s 

2Y  P 
o  o 


(10) 


This  expression  is  plotted  in  Fig.  4,  using  the  Ys  of  Eq.  (9). 
For  an  ideal  gas  of  Y  =  1.4,  it  reduces  to 


U 

s 


V 


AP 


S  ( Y  +  1  ) 
2Y 


1.087  /0.0583AP  +  1  kft/s  ,  (11) 

s 


in  which  C0  =  /YP0/p0  =  1.087  kft/s  and  P0  for  sea  level  =14.7  psi. 

In  air  at  sea  level,  this  simple  form  [Eq.  (11)]  is  within  5  percent 
correct  for  APS  <  1,000,000  psi.  As  can  be  seen  in  the  comparison 
plot  of  Fig.  4,  the  actual  shock  velocity  is  only  slightly  lower  than 
the  ideal  gas  simple  form,  and  the  dependence  of  shock  velocity  on  the 
equation  of  state  for  air  is  minimal.  For  most  purposes,  this  simple 
form  is  sufficiently  accurate,  since,  in  reality,  the  local  variations 
in  fireball  growth  far  exceed  5  percent.  Above  =200,000  psi,  the 
front  velocity  is  increased  by  radiative  effects,  and,  therefore,  the 
Hugoniot  expression  is  inadequate. 
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Figure  4.  Shock  velocity  versus  shock  overpressure  in  air  at  sea  level. 


PEAK  PARTICLE  VELOCITY. 

The  air  speed  immediately  behind  an  air  shock  can  be  written  in 
terms  of  shock  overpressure  [from  Eqs.  (1),  (2),  and  (3)]  as 


u 

s 


2AP 

_ s 

Y  P 
o  o 


AP  +  P 
s  o 


DAP  +  2Y  P  ] 
s  so 


1/2 


(12) 


For  an  ideal  gas  of  Ys  =  YQ 


u 

s 


C  AP 
o  s 


7yp  [ (Y  +  DAP  + 
o  s 


=  1.4,  this  expression  simplifies 

0 . 21 88 AP 

-  =  - - —  kft/s 

2YP  ] / 2  JT?  +  17.15 

o  s 


to 


(13) 


This  ideal  gas  form  is  valid  for  shocks  in  air  to  within  5  percent  for 
APS  <  1,000,000  psi.  The  actual  particle  velocity  is  higher  by  about 
4  to  5  percent  between  10,000  and  50,000  psi,  and  by  less  than  that 
elsewhere. 

For  strong  shocks,  i.e.,  for  APS  >>  P0, 

u  =  C  /2AP  / [ P  Y  ( Y  +  1 ) ]  •  (14) 

s  o  s  o  o  s 

[Equation  (9)  can  be  used  to  find  Ys  as  a  function  of  APS.] 

The  peak  particle  velocity  as  defined  by  Eqs.  (12)  and  (9)  is 
plotted  in  Fig.  5,  and  compared  to  the  ideal  gas  form  of  Eq.  (13). 

PEAK  DYNAMIC  PRESSURE. 

The  blast  wind  pressure  or  dynamic  pressure  Q  is  defined  as 


At  the  shock  front,  using  Eqs.  (7)  and  (12),  the  peak  dynamic 
pressure  Qs  becomes 
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Q„  = 


For  an  ideal  gas  (Ys  =  Y0  =  1.4)  and  standard  conditions, 


Q 


s 


2YP 

o 


+  (Y  -  1 )AP 

s 


41.2+  0. 4AP 

s 


psi  . 


(17) 


A  better  approximation  for  real  air  is 


ws  =  tt(  1  +  0.241  tt  +  0.437 6tt  ) _ 

AP  2  7 

s  0.041  +  0 . 4tt  +  0 . 02891  tt  +  0.1015tt 


0.01  251  it 


7.29x1  0  8tt4 


9.763x1  0  10tt 


1  +  9.649x1  0  7tt5  1  +  2.61x10  2V2  1  +  6.957x1  0  28tt12 


—  ft  ft 

5.052x10  tt 


5  6 

6.021 xl  Q  c 


,  12  ,  —  r~ 1 1 1  J  2  1  2 

1  +  1.368x1  0  tt  1  +  3.541x10  c 


9  17  in8r14 

2.17x10  C 


0.7670c 


2.839 


1  +  1.62xl09c15  1  +  0.1646c8*678 


(18) 


where  tt  =  APS/1000  (APS  in  pounds  per  square  inch), 
C  =  tt/1000. 


The  above  expression  is  accurate  to  better  than  3  percent  for  all 
overpressures  less  than  7,000,000  psi.  If  the  last  term  is  replaced 
by 


0.6805c8 

4  ’ 

1  +  0.1230c 


(19) 
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the  expression  is  accurate  to  better  than  4  percent  for  overpressures 
less  than  6,500,000  psi.  The  standard  of  comparison  is  Eq.  (16)  using 
Eq.  (9)  for  Ys.  Again,  above  about  200,000  psi,  for  nuclear  blast 
waves,  the  influences  of  radiation  transport  are  likely  to  invalidate 
this  adiabatic  shock  expression  for  peak  dynamic  pressure. 

The  ratio  of  peak  dynamic  pressure  to  peak  overpressure  is 
plotted  in  Fig.  6,  and  compared  with  the  fit  from  Eq.  (18).  The 
figure  shows  that  the  quadratic  dependence  of  dynamic  pressure  on 
overpressure  at  low  overpressures  gives  way  to  something  nearer  linear 
dependence  above  100  psi.  Shock  dynamic  pressure  is  equal  to  shock 
overpressure  at  around  70  psi. 

SHOCK  TEMPERATURE. 

The  increased  temperature  in  a  shock  front  A9S  can  be  expressed, 
using  the  gas  law  relation  P  =  p$0  and  the  Hugoniot  equations 
[Eqs.  (1),  (2),  and  (3)]  as: 


in  which  $0  is  the  gas  constant  for  air  at  standard  conditions,  $s  is 
the  gas  constant  for  shocked  air ,  and  9q  is  the  ambient  (preshock) 
temperature.  For  an  ideal  gas,  Y  =  Y0  =  Ys,  and  10  =  $g.  The  tem¬ 
perature  ratio  becomes: 


(21  ) 


For  standard  conditions  in  air,  Y 
and 


1.4, 


PQ  =  14.7  psi,  90 


273-2  K, 
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versus  shock  overpressure. 


5.31 AP  (1  +  0.034AP  ) 

_ s _ s__ 

( 1  +  0.0583AP  ) 

3 


°c  . 


(22) 


A  better  approximation  to  the  shock  temperature  increase  in  air 
is  given  by: 


A  0 


531  Ott(  1  +  34tt) 


1  +  58. 3 tt  + 


6. 53tt 


3.5 


(1  +  0.2027tt  ) 


'C. 


(23) 


A  still  better  approximation  to  the  temperature  of  a  shock  (in 
air  at  sea  level,  standard  conditions)  combines  Eq.  (20)  with  the 

following  fit  to  the  thermal  equation  of  state  on  the  Hugoniot: 

) 


0. 0081  68-tt2  +  0. 00192*1  tt2 

(1  +  0.01807tt2)  (  1  +  0.0005596tt2) 


+ 


78.45c2  +  0.0009018c14 

1  +  10.12c2  1  +  0.0001 484c4 


(24) 


where  tt  =  APS/1000, 

C  =  it  / 1  000  =  APsxl0-6. 

THERMAL  EQUATION  OF  STATE  FOR  AIR. 

Actually,  the  gas  "constant"  $s  increases  manyfold  as  the  shock 
overpressure  rises,  so  that  at  high  pressures  the  real  air  shock 
temperature  falls  below  the  ideal  gas  temperature  by  factors  ranging 
up  to  more  than  16. 

A  relation  between  the  gas  constant  %  and  temperature  and  density 
allows  the  evaluation  of  Eq.  (20)  and  the  determination  of  the  shock 
temperature  as  a  function  of  the  shock  overpressure. 

A  fit  to  the  thermal  equation  of  state  for  air,  based  on 
properties  of  air  derived  by  Gilmore  [1955,  1959]  and  Hilsenrath, 
Green,  and  Beckett  [1957]  is  as  follows: 
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0.8dS 

1  +  dB 
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%_  pv_ 

%0  =  Xoe 


4aa 


8 


1  +  aa 


8 


8. 4ba  |  2ca 

o  8 

1  +  ba  1  +  ca 


0.2f  B 
+  fB  ’ 

(25) 


where  a  -  10  4x0q  0*086  for  q  degrees  kelvin,  n  =  p/p0, 


B 

a 

b 

c 
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a 


0.9746  +  0.0254n~°*21 556 
0.7778x1 o-1 6, 

0.602x1 o~8  t 

0.5097 xl 0~8  f 


f 


d  =  2.20, 
f  =  0.971 xio4, 


with  =  0.04161  for  pressure  in  pounds  per  square  inch,  air  tempera¬ 
ture  6  in  degrees  kelvin,  density  p  in  kilograms  per  cubic  meter,  and 
specific  volume  V  =  1/p. 

The  thermal  equation  of  state  is  illustrated  in  Fig.  7.  The 
figure  shows  a  plot  of  PV/$06  versus  temperature  0,  which  is  the 
equivalent  of  the  ratio  of  the  gas  constant  to  that  for  ambient  sea- 
level  ai r  %/%0.  It  is  also  a  measure  of  the  number  of  free  particles 
per  mole  of  air  at  a  temperature  0  and  a  density  ratio  n  relative  to 
the  number  in  a  mole  of  air  at  standard  sea-level  conditions.  At  a 
little  less  than  10,000  K,  the  oxygen  and  nitrogen  molecules  in  air 
dissociate,  and  the  relative  number  of  particles  rises  by  a  factor  of 
two.  At  extreme  temperatures,  like  1  million  degrees,  all  the 
electrons  are  freed  from  the  atoms,  and  the  average  number  of  par¬ 
ticles  (ions  and  electrons)  rises  to  more  than  16  per  original  air 
molecule. 

Again,  the  changes  in  the  gas  constant  ratio  v  are  slow  relative 
to  density  and  pressure  (or  temperature)  changes,  so  that  simple 
iterative  methods  can  be  used  to  determine  the  proper  thermodynamic 
state. 
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Air  temperature  (K) 

Figure  7.  Thermal  equation  of  state  for  air  versus  temperature  and  density. 


NORMAL  REFLECTION  OF  SHOCKS. 

When  a  shock  wave  strikes  a  rigid  plane  head  on  (as  illustrated 
in  Fig.  8),  the  reflected  shock  conditions  can  be  derived  from  conser¬ 
vation  considerations  similar  to  those  expressed  in  Eqs.  (1),  (2),  and 
(4),  for  mass,  momentum,  and  energy  across  the  reflected  shock: 


PRUR  =  (U5  +  UR)PS 


PR  -  Ps  '  (us  *  UR)2pa  -  UR2pR  ’  UsURpR 


R 


p  +  p  / 

R  s  /  1 


Or^r"15 


Ps  PR 


(mass) ,  (26) 

(momentum),  (27) 


(energy) ,  (28 ) 


in  which  subscripts  R  refer  to  conditions  after  reflection  and  sub¬ 
scripts  s  apply  to  shock  values  prior  to  reflection. 
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7 
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Figure  8.  Reflected  shock  conditions. 
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Equations  (26),  (27),  and  (28)  together  with  the  usual 
tions  [Eqs.  (1),  (2),  and  (4)],  lead  to  a  quadratic  in 
factor  RF  =  (APR/APS): 


Hugoniot  rela¬ 
the  reflection 


[ 2Y  P  + 
R  o 


[2Y  P  + 
0  0 


-  1 )AP  ] 

-  1)APS] 


=  0  . 


(29) 


For  an  ideal  gas  (where  YR  =  Ys  =  Y0  =  Y),  that  expression  simplifies 
to: 


RF 


AP 

s 


2  +  (Y  +  1  ) 


(30) 


in  which  APR  is  the  reflected  peak  overpressure  from  a  normally  inci¬ 
dent  shock  of  peak  overpressure  APS  in  an  ambient  atmosphere  of  pres¬ 
sure  P0,  and  Qs  is  given  by  Eqs.  (16),  (17),  or  (18).  Equation  (18), 
though,  is  an  expression  for  Q^/ APg  for  a  shock  in  real  (sea— level) 
air,  and  is  therefore  inconsistent  with  the  ideal  gas  assumption  of 
Eq.  (30). 

For  air  at  sea  level,  a  better  approximation  to  the  reflection 
factor  RF  is  given  by  the  formula 
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1 


RF  =  2  + 


2 . 6  5  5tt 


1  +  0.  1  7  2  8  tt  +  0.001921 


TT 


+  0.004218  4-  48 . 34tt  +  6. 856-tt' 

1  +  7  - 9 9 7 tt  +  3  .844-rr2 


(31  ) 


where  tt  -  APS/1000,  (APS  in  pounds  per  square  inch). 


For  a  strong  shock  ( APS  >>  P0),  the  expression  in  Eq.  (29) 
simplifies  to 


2TR 

RF  =  1  +  - — ~  .  (32) 

s 

In  Fig.  9,  the  reflection  factor  predicted  by  the  approximation 
in  Eq.  (31)  is  compared  with  a  more  exact  solution  from  Eq.  (29)  using 
the  equation  of  state  for  air  given  in  Sec.  5  [Eqs.  (69)  through 
(74)]. 

TABLE  OF  SHOCK  VARIABLES. 

The  shock  parameters  treated  in  this  section  (and  plotted  in 
Figs.  2  through  7,  and  9)  are  listed  for  a  range  of  shock  overpres¬ 
sures  in  Table  1,  for  standard,  sea-level  air.  The  specific  heat 
ratio  Ys  is  from  Eq.  (9),  the  shock  density  ratio  ns  from  Eqs.  (7)  and 
(9),  the  shock  velocity  Us  from  Eqs.  (9)  and  (10),  the  peak  particle 
velocity  us  from  Eqs.  (9)  and  (12),  and  the  peak  dynamic  pressure  Qs 
from  Eqs.  (9)  and  ( 1 6 ) .  The  ratio  of  the  gas  constants  (thermal)  v  is 
given  by  Eq.  (25);  the  shock  temperature  rise  A6S  is  provided  by 
combining  Eqs.  (9),  (20),  and  (25),  and  the  normal  reflection  factor 
RF  is  arrived  at  from  Eq.  (29)  and  the  equat ion-of-state  for  air,  from 
Sec.  5,  Eqs.  (69)  through  (74). 
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Figure  9.  Normal  reflection  factor  versus  shock  overpressure  for  sea-level  air. 


Table  1.  Shock  variables. 


APS  ,  _  ,  Is  us  us  Q 


(psi) 

'S 

(PS/P0) 

(kft/s) 

(kft/s) 

(psi) 

5 

(°C) 

(APr/APs) 

.1 

.400 

1.005 

1.090 

.00527 

.000242 

1.00 

0.5 

2.009 

.15 

.4 

1.007 

1.092 

.00789 

.000545 

1.00 

0.8 

2.012 

.2 

.4 

1.010 

1.094 

.01051 

.00097 

1.00 

1.1 

2.014 

.3 

.4 

1.015 

1.097 

.01571 

.00218 

1.00 

1.6 

2.019 

.4 

.4 

1.019 

1.100 

.02089 

.00387 

1.00 

2 

2.025 

.6 

.4 

1.029 

1.106 

.0312 

.00869 

1.00 

3 

2.035 

.8 

.4 

1.039 

1.112 

.0413 

.01541 

1.00 

4 

2.045 

1 

.4  . 

1.048 

1.118 

.0514 

.0241 

1.00 

5 

2.055 

1.5 

.4 

1.072 

1.134 

.0760 

.0539 

1.00 

3 

2.080 

2 

.4 

1.095 

1.149 

.1000 

.0953 

1.00 

10 

2.105 

3 

.4 

1.146 

1.178 

.1462 

.2125 

1.00 

15 

2.154 

4 

.4 

1.187 

1.207 

.1903 

.3742 

1.00 

20 

2.202 

6 

.4 

1.275 

1.263 

.2728 

.8264 

1.00 

28 

2.297 

7 

.4 

1.318 

1.290 

.  3116 

1.1146 

1.00 

33 

2.343 

10 

.4 

1.443 

1.368 

.4200 

2.214 

1.00 

45 

2.478 

15 

.4 

1.636 

1.489 

.5788 

4.771 

1.00 

64 

2.692 

20 

.4 

1.814 

1.600 

.7179 

3. 137 

1.00 

82 

2.891 

30 

.4 

2.129 

1.803 

.9559 

16.93 

1.00 

117 

3.254 

40 

.4 

2.400 

1.985 

1.158 

27.99 

1.00 

150 

3.575 

50 

.399 

2.640 

2.150 

1.336 

41.01 

1.00 

132 

3.362 

70 

.397 

3.043 

2.448 

1.643 

71.50 

1.00 

244 

4.351 

100 

.392 

3.522 

2.833 

2.028 

126.1 

1.00 

332 

4.930 

150 

.385 

4.092 

3.378 

2.551 

231.9 

1.00 

475 

5.629 

200 

.378 

4.508 

3.844 

2.990 

350.8 

1.00 

612 

6.126 

300 

.368 

5.063 

4.637 

3.718 

609.4 

1.00 

882 

6.795 

450 

.355 

5.602 

5.612 

4.608 

1035.5 

1.00 

1268 

7.413 

700 

.336 

6.207 

6.927 

5.807 

1322 

1.00 

1866 

3.043 

1000 

.312 

6.831 

8.207 

7.002 

2916 

1.00 

2486 

8.559 

1500 

.279 

7.728 

9.953 

3.660 

5046 

1.00 

3359 

9.215 

2000 

.252 

8.569 

11.409 

10.074 

7569 

1.01 

4045 

9.758 

3000 

.237 

9.175 

13.91 

12.39 

12263 

1.10 

5273 

10.656 

4000 

.228 

9.565 

16.02 

14.34 

17130 

1.18 

6313 

11.375 

6000 

.207 

10.51 

19.52 

17.66 

28530 

1.24 

3288 

12.430 

3000 

.191 

11.35 

22.45 

20.47 

41390 

1.34 

9530 

13.132 

10000 

.180 

11.99 

25.04 

22.94 

54940 

1.46 

10340 

13.604 

15000 

.172 

12.56 

30.60 

28.16 

86660 

1.31 

12030 

14.20 

20000 

.179 

12.12 

35.39 

32.46 

111200 

2.08 

14450 

14.35 

30000 

.207 

10.65 

43.61 

39.51 

144700 

2.62 

19750 

14.13 

40000 

.210 

10.50 

50.39 

45.59 

190000 

3.12 

22400 

13.66 

60000 

.218 

10.16 

61.83 

55.74 

274800 

4.04 

26900 

12.64 

80000 

.234 

9.537 

71.65 

64.13 

341500 

4.76 

32500 

11.  72 

100,000 

.232 

9.612 

80.07 

71.74 

430600 

5.19 

37000 

10.95 

150,000 

.209 

10.56 

97.55 

88.32 

717200 

5.90 

44500 

9.520 

200,000 

.212 

10.43 

112.7 

101.9 

942900 

6.74 

52600 

3.552 

300,000 

.255 

8.841 

139.4 

123.5 

1176000 

8.73 

72000 

7.341 

400,000 

.2545 

8.857 

160.9 

142.8 

1571000 

9.75 

85800 

6.619 

600,000 

.254 

8.873 

197.1 

174.9 

2362000 

11.0 

113700 

5.800 

300,000 

.264 

8.575 

228.1 

201.5 

3030000 

11. 7 

148000 

5.346 

1000,000 

.285 

8.017 

256.2 

224.2 

3508000 

12.1 

192000 

5.060 

1500,000 

.353 

6.665 

318.4 

270.0 

4249000 

12.3 

339000 

4.659 

2000,000 

.400 

6.000 

371.3 

309.4 

5000000 

12.4 

500000 

4.449 

3000,000 

.391 
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SECTION  3 

FREE-AIR  AND  SURFACE  NUCLEAR  BURSTS 


PEAK  OVERPRESSURE  VERSUS  RANGE  (FREE-AIR  BURST). 

An  expression  that  is  appropriate  for  overpr essures  in  the  range 
0.07  <  APS  <  400,000  psi  is 


AP 


s 


2  3 

—  +  - - —  + 

r  1  .5 
r 


1 05 ( 1  Or) 


1  +  1  2 , 920( 1  Or) 


8 


ps  1 


(33) 


where  r  is  scaled  range  in  kilofeet  per  cube-root  kiloton.  In 
general,  the  agreement  with  a  detailed  one-dimensional  calculation 
[Brode,  1959b]  is  within  5  percent.  Extension  to  pressures  below 
2  psi  is  in  agreement  with  the  "Combined  Airborne  Polynomial"  to 
within  5  percent,  as  well.  The  term  in  brackets  is  unnecessary  for 
overpressures  below  200  psi  (viz.,  r  >  0.22  kft/KT1^).  This  fit  is 
only  10  percent  high  at  600,000  psi.  At  the  lower  extreme,  it  is 
20  percent  high  at  0.01  psi  in  comparison  with  the  suggested  curve, 
which  is  an  extension  of  the  Airborne  curve  and  other  (theoretical) 
extrapolations.  The  fit  is  in  agreement  with  the  1 -KT  standard  [Need¬ 
ham  and  Crepeau,  1981]  to  within  10  percent  (beyond  20-m  scaled 
range).  The  fit  is  plotted  in  Fig.  10. 

The  accuracy  of  the  fit  should  not  mislead  the  user  into  believ¬ 
ing  that  such  overpressures  are  repeatable  or  reproducible  with  com¬ 
parable  accuracy.  A  review  of  the  original  atmospheric  nuclear  test 
measurements  suggests  a  much  larger  data  scatter  even  for  the  control¬ 
led  test-site  conditions  and  surfaces.  The  scatter  in  test  data  far 
exceeds  the  few  percent  difference  between  "best-guess"  curves  and 
fits.  Figure  1 1  is  a  graph  of  approximate  data  scatter,  as  derived 
from  Brode  [1981],  plotted  against  peak  overpressure.  At  the  higher 
overpressure,  the  plot  suggests  a  minimum  uncertainty  of  70  percent 
(from  minimum  to  maximum  datum) ,  which  translates  into  a  pressure 
variation  of  a  factor  of  5.  The  uppermost  curve  in  Fig.  11  roughly 
corresponds  to  a  21  limit;  i.e.,  bounding  about  95  percent  of  the 
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Peak  overpressure  (psi) 


Figure  10.  Peak  overpressure  versus  shock  radius  for  1-KT 

free-air  nuclear  burst  in  standard,  sea-level  air. 
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Figure  11.  Scatter  in  percent  of  range  to  given  peak  overpressure 
for  near-surface  atmospheric  nuclear  tests. 


data.  The  lower  curve  approximates  a  IE  value.  The  circles  are 
estimates  of  the  IE  value  from  the  plots  in  Brode  [1981],  while  the 
arrows  represent  a  value  of  approximately  2E. 

PEAK  OVERPRESSURE  VERSUS  RANGE  (IDEAL-SURFACE  BURST). 

The  free-air  curve  can  be  recast  for  a  rigid  reflecting  surface 
by  using  2W  in  place  of  yield  W  in  the  scaled  ground  range  x: 


AP 


s 


2.52  +  4.2^ 

x  1  .5 

x 


+  41  .7Q0X)11 

1  +  2035(  1  Ox) 


psi, 


(34) 


where  x  is  in  kilofeet  per  cube-root  kiloton. 


TEST-SURFACE  PEAK  OVERPRESSURE  VERSUS  RANGE. 

Overpressures  from  near-surface  bursts  in  early  atmospheric 
nuclear  tests  showed  peaks  generally  below  those  of  ideal  surface 
bursts  that  can  be  approximated  by  scaling  1 .8W  rather  than  2W.  The 
observed  lower  peaks  were  below  1000  psi  and  mostly  above  8  ps i .  In 
that  range,  the  approximation  (using  1 .8W  scaling)  becomes 

AP  »  —  +  4—  «■  2-A-  psi,  for  6  <  AP  <  1000  psi.  (35) 

s  x  x1.5  x3  5 


SHOCK  RADIUS  VERSUS  PEAK  OVERPRESSURE  (FREE-AIR  BURST). 

The  following  fit  can  be  used  to  search  for  the  range  to  a  given 
shock  overpressure.  Employing  this  fit  saves  one  from  the  tedium  of 
inverting  the  approximation  of  Eq.  (33)  for  overpressure  as  a  function 
of  range.  The  fit  is  accurate  to  within  2  percent  of  Eq.  (33): 


sr 


m 


2.463  + 

0.9846 


1  .176 


AP 


1/3 


0 . 0004726tt 


1  .5 


—  /I  q 

1  +  2.952x1 0~ 


kft 


(36) 


This  approximation  is  for  a  free-air  burst  in  which  sr  is  the  range  in 
kilofeet;  m,  the  scaling  factor,  is  the  cube-root  yield  in  cube-root 
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kilotons;  APS  is  in  pounds  per  square  inch;  and  it  =  APS/1000.  (Note 
that  sr  is  equal  to  m  *  r.)  The  comparable  range  for  an  ideal-surface 
burst  ( 2W)  is  2^3  (-1.26  times  larger).  For  bursts  over  the  (empiri¬ 
cal)  test-site  surface  (1.8W),  the  range  is  -1.22  times  larger  than 
that  for  a  free-air  burst. 

PEAK  OVERPRESSURE  VERSUS  TIME  OF  ARRIVAL  (FREE-AIR  BURST). 

The  peak  overpressure  can  also  be  expressed  as  a  function  of  the 
scaled  time  of  arrival: 


psi, 


(37) 


where  i  5  T/m, 

m  s  (in  cube-root  kilotons), 


T  =  time  of  arrival  (in  milliseconds). 

This  fit  to  the  time  of  arrival  and  peak  overpressure  is  based  on 
early  detailed  calculations  [Brode,  1959b,  1966]  scaled  to  1  KT  for 
pressures  between  2  £  APS  £  500,000  psi.  The  fit  is  good  to  better 
than  4  percent  below  50,000  psi,  and  to  better  than  3  percent  below 
10,000  psi. 

PEAK  OVERPRESSURE  VERSUS  TIME  OF  ARRIVAL  (SURFACE  BURST). 

Using  the  2W  assumption  for  surface  bursts,  the  fit  from 
Eq.  (37)  becomes 


916,000  +  570t 


0.75 


+ 


psi. 


(38) 
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TIME  OF  ARRIVAL  VERSUS  PEAK  OVERPRESSURE. 

Conversely,  the  time  of  arrival  T  can  be  expressed  as  a  function 
of  the  peak  overpressure  (for  a  free-air  burst)  as: 

T  ~  m  |o. 03394  + 

s 

where  m  =  W1 ^3  (in  cube-root  kilotons), 

W  s  yield  (in  kilotons), 

APS  =  pounds  per  square  inch. 

For  a  surface  burst,  one  uses  2W.  The  above  form  is  accurate  to 
±2  percent  for  2  ^  APS  £  100,000  psi.  This  relation  for  a  1-KT,  free- 
air  burst  is  illustrated  in  Fig.  12. 

As  is  evident  in  Fig.  12,  the  time  of  arrival  is  nearly  inversely 
proportional  to  the  overpressure.  Actually,  the  product  of  time  of 
arrival  and  peak  overpressure  to  the  fractional  power  0.875  is  a 
slowly  varying  function,  and  can  generate  a  more  readable  curve.  That 
relation  is  plotted  in  Fig.  13.  The  approximate  proport ionali ty  means 
that  the  product  of  overpressure  times  arrival  time  to  the  1.14  power 
is  roughly  constant  (APS  x  T1*^  ^constant). 

TIME  OF  ARRIVAL  VERSUS  SHOCK  RADIUS. 

The  time  of  arrival  T  can  also  be  represented  in  terms  of  the 
shock  radius  sr  from  a  free-air  burst  [Brode,  1970], 

_  0.54291m3  -  21.l85(sr)m2  +  36l.8(sr)2m  +  2383(sr)3  ^  (4o) 

m2  +  2. 048( sr )m  +  2.6872(sr)2 

with  sr  in  kilofeet,  m  =  W1/3  (in  cube-root  kilotons),  and  W  =  yield 
( in  kilotons) . 

This  shock  arrival-time  form  is  limited  to  the  range  of  the 
"Empirical  59"  data  [Moulton,  I960],  namely,  620  ys  to  26  s  at  1  KT, 


893 


AP 


0.80424 


2015 

AP 


ms , 


(39) 
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Scaled  time  of  arrival  (ms/KT 


Peak  overpressure  (psi) 
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Peak  overpressure  (psi) 


Figure  12.  Time  of  arrival  versus  peak  overpressure  for  1-KT 
free-air  burst  at  sea  level. 
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Peak  overpressure  (psi) 

Product  of  time  of  arrival  and  peak  overpressure^’**^  versus 
peak  overpressure  for  1  KT  (T  in  milliseconds). 


which  spans  a  range  of  free-air  peak  overpressures  from  17,200  to 
0.07  psi.  (Use  2W  in  m  for  a  surface  burst;  i.e.,  replace  m  with 
m  x  2^3 . ) 


CLOSE-IN  TIME  OF  ARRIVAL  VERSUS  SHOCK  RADIUS. 

For  higher  pressure  (or  earlier  time  of  arrival  T),  the  following 
fit  is  appropriate  [Brode,  1970]: 

T  =  — —  ms,  (41  ) 

b  +  c 

where  a  =  [0.543m2  -  21.8(sr)m2  +  386(sr)2m  +  2 383 ( sr ) 2]sr8 , 
b  =  [2.99*10~8m6  -  1  .91  *1  0“4(  sr  )2m4  +  1.032(sr)4m2 
-  4. 43( sr )8]m41 0“8, 

c  =  [1.028m2  +  2.087(sr)m  +  2.69(sr )2]sr8 . 

Equation  (41)  is  alternatively  written  in  terms  of  the  scaled  range  r 
in  kilofeet  per  cube-root  kiloton: 


where  a  =  (0.543  -  21 .8r  +  386r2  +  2383r2)r8, 

b  =  (2.99x10*8  -  1  .91x1 0“V2  +  1.032r4  -  4.43r6)10~6, 
c  =  (1.028  +  2. 087r  +  2.69r2)r8, 

and  is  valid  for  10“3  <  T  <  26,000  ms  at  1  KT.  This  more  complex  fit 
is  advisable  for  pressures  above  10,000  psi,  or  scaled  times  less  than 
0.6  ms/KT1^.  This  relation  is  plotted  in  Fig.  14. 

For  a  surface  burst,  use  m  =  (2W)^3. 

POSITIVE  OVERPRESSURE  DURATION  VERSUS  TIME  OF  ARRIVAL  (FREE-AIR  BURST). 

The  duration  of  the  positive  phase  for  overpressure  can  be  ex¬ 
pressed  variously  as  a  function  of  overpressure,  shock  range,  or  time 
of  arrival  in  a  fit  to  detailed  one-dimensional  blast  calculations 
[Brode,  1959b]: 
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Scaled  time  of  arrival  (ms/KT l/J) 


Scaled  shock  radius  (ft/KT^) 


Figure  14. 


Time  of  arrival  versus  shock  radius  for  1-KT  free-air 
burst. 
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Scaled  time  of  arrival  (ms/KT 


D 

P 


(m)  (81 3,1 40m2  +  1 1 ,41 2mT  +  31 3T2) 
(6,780m2  +  444. 7mT  +  T2) 


(43) 


where  T  is  in  milliseconds. 

POSITIVE  OVERPRESSURE  DURATION  VERSUS  TIME  OF  ARRIVAL  (SURFACE  BURST). 
For  a  surface  burst,  the  form  in  Eq.  (43)  becomes 

D+  -  («)  LU62J,30.0.n^.J.8,115mT2  »  3W2),  ms>  W) 
p  (10,790m2  +  560. 3mT  +  T2) 

where  T  is  in  milliseconds.  It  is  assumed  that  the  effect  of  the 
surface  is  to  double  the  blast  yield. 

POSITIVE  OVERPRESSURE  DURATION  VERSUS  PEAK  OVERPRESSURE  (FREE-AIR 
BURST) . 

Expressed  as  a  function  of  peak  overpressure ,  the  duration  of 
positive  overpressure  can  be  approximated  by 


D 

P 


(m) 


-148.6  + 


497.3 


1  629ir 


0.871 1 


(1 


18.68tt0,6783) 


(1  +  6. 477tt°*8555) 


ms , 

(45) 


where  it  =  APS/1000,  in  kilopounds  per  square  inch. 

This  form,  in  combination  with  Eq.  (33),  is  compared  in  Fig.  15  with 
both  the  calculation  to  which  it  is  a  fit  [Brode,  1959b]  and  the  DNA 
1-KT  standard  values  [Needham  and  Crepeau,  1981].  The  difference 
between  those  two  predictions  is  typical  of  the  uncertainty  and  ex¬ 
pected  variation  (particularly  at  large  distances)  in  blast  parameters 
such  as  positive  phase  durations.  For  a  surface  burst,  this  form 
should  be  multiplied  by  2173  =  1.26  (i.e.,  use  2W  in  m). 
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Figure  15.  Fit  compared  to  calculations:  overpressure  durations  at  1  KT. 


POSITIVE  OVERPRESSURE  DURATION  VERSUS  RANGE. 

As  a  function  of  scaled  range  (1  KT,  free-air  burst),  the  dura¬ 
tion  of  positive  overpressure  can  be  approximated  by- 

ms, 

(46) 


D  =  (m) 
P 


69.12  + 


46.19 


4043r 


6.329 


(1  +  3,000,000r7’217) 


(1  +  37.l6r5*621  ) 


where  r  is  scaled  range  in  kilofeet  per  cube-root  kiloton.  Again,  for 
a  surface  burst,  Dp  should  be  increased  by  a  factor  of  21/3,  i.e.,  the 
free-air  form  for  twice  the  yield. 


OVERPRESSURE  IMPULSE  IN  POSITIVE  PHASE  VERSUS  PEAK  OVERPRESSURE  (FREE- 
AIR  AND  SURFACE  BURST). 

Positive  phase  impulse  Ip  is  defined  by  the  integral 


I 

P 


T+D 

P 

AP(t)  dt 
T 


m 


t  +  D 

/ 


T 


AP(o)  da  . 


(47) 


A  simple  approximation  to  this  impulse  as  a  function  of  peak  overpres¬ 
sure  is 


I  *  "  TTp-  m  psi-ms.  (48) 

H  (1  +  0.00385AP  ) 

s 

This  fit,  when  used  with  Eq.  (33) »  leads  to  the  values  compared 
in  Fig.  16  versus  radius.  The  impulse  approaches  a  constant  at  small 
ranges  and  decays  approximately  as  the  inverse  of  the  range  elsewhere. 
For  a  surface  burst,  that  expression  should  be  multiplied  by  1.26 
(i.e.,  by  2^3),  which  leads  to  replacing  the  coefficient  145  by  183. 
This  form  is  good  to  better  than  1 0  percent  for  2  <  APQ  <  1 00  000 
psi.  Comparison  between  the  approximation  as  a  function  of  peak 
overpressure  and  the  detailed  numerical  results  [Brode,  1964]  is  made 
in  Fig.  17.  In  that  plot,  it  is  evident  that  impulse  increases 
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Fit  compared  to  detailed  calculation  and  1-KT  standard:  over¬ 
pressure  impulse  in  positive  phase  (1-KT,  free-air  burst). 
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Figure  17.  Fit  compared  to  calculation:  impulse  in  positive  overpressure 
versus  peak  overpressure  (1-KT,  free-air  burst). 


approximately  as  the  cube  root  of  overpressure ,  but  tends  toward  a 
constant  at  about  a  few  thousand  pounds  per  square  inch. 

OVERPRESSURE  VERSUS  TIME. 

The  following  analytic  expression  is  valid  for  overpressures  less 
than  about  15,000  psi.  It  is  an  approximate  form,  modified  from 
earlier  fits  [Brode,  1970,  1978]  for  the  overpressure  in  the  positive 
phase  as  a  function  of  time.  In  these  fits,  time  is  zero  at  the 
instant  of  burst: 


AP  = 


AP  (a) 
s 


I 


0.417  +  0.583 


40  +  t2 


c-t\  . 

1  /  P31’ 

D  /m  / 

P 

(49) 


where  APs(o)  is  found  in  Eq.  (37)  and  o  replaces  t,  t  =  T/m,  a  =  t/m, 
T  is  the  time  of  arrival  (in  milliseconds),  t  is  the  time  after  time 
of  arrival  (t  Z  T,  i.e.,  a  £  t) .  Both  T  and  t  are  measured  from  the 
instant  of  detonation.  The  scale  factor  m  =  (in  cube-root 

kilotons)  and  Dp  is  the  duration  of  the  positive  phase  [Eqs.  (43) 
through  (46)].  As  before,  for  a  surface  burst,  use  m  =  (2W)1/3, 


OVERPRESSURE  VERSUS  TIME  (SURFACE  BURST). 

Alternatively,  one  can  use  the  zero  burst  height  from  the  more 
complex  equations  for  height  of  burst  and  range  [Eq.  (63)  and  Speicher 
and  Brode,  1980a, b,  1981,  1984a, b].  When  zero  burst  height  is  in¬ 
serted  and  expressions  are  simplified,  the  calculation  for  a  surface 
burst  becomes 


AP  =  AP  (t) 
s 


(50) 
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where 


D  = 
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-9  S  -10  4  2R 

1  +  1.473x1  0  x  1  +  3.207x1  0  iuth*^ 


0.03209X 


1.25 


0. 1 966x 


1 .22 


1  +  9.914x10”8x4  1  +  0.767x1,22 


and  X  is  the  scaled  ground  range  in  feet  per  cube-root  kiloton. 

Note  that  in  this  form,  unlike  Eq.  (49),  the  peak  overpressure  is 
the  peak  at  the  time  of  shock  arrival.  The  peak  overpressure  (APS) 
for  this  surface  burst  form  can  be  approximated  in  various  ways. 

Using  Eq .  (38),  which  is  Eq.  (37)  with  2W,  leads  to  a  convenient  form 
in  terms  of  arrival  time. 

Using  y  =  0  in  the  height  of  burst  form  from  Speicher  and  Brode 
[1980b,  as  revised  in  1981  and  further  revised  in  1984a]  (i.e.,  for  a 
surface  burst)  leads  to  a  function  of  scaled  ground  range  somewhat 
different  than  the  form  of  Eq.  (34),  namely: 
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4.166 


-  0.2905 


Ap  =  10-4?  +  2.9902 
s  -  X1 .22  x3.053 


+ 


(51  ) 


where  x  is  scaled  ground  range  in  kilofeet  per  cube-root  kiloton. 

Equation  (34)  is  preferred  over  Eq.  (51).  Eq.  (34)  is  a  bit 
simpler  algebraically  and  is  also  a  closer  fit  to  the  detailed  one¬ 
dimensional  calculation  [Brode,  1959b].  Equation  (51)  is  within 
7  percent  of  Eq.  (34)  between  10,000  psi  and  0.03  psi.  The  largest 
deviation  occurs  near  150,000  psi,  a  point  in  the  detailed  calculation 
just  after  the  air  shock  has  fully  formed  and  the  debris  shock  has 
overtaken  it  [at  that  point,  Eq.  (51)  is  low  by  -23  percent).  At 
closer  ranges  and  higher  pressures,  details  of  the  nuclear  weapon 
become  increasingly  important,  so  that  no  general  curve  fit  will  be 
wholly  correct.  In  addition,  radiation  transport  plays  an  increas¬ 
ingly  important  role,  making  normal  shock  physics  less  applicable. 

DURATION  OF  OUTWARD  DYNAMIC  PRESSURE  VERSUS  PEAK  OVERPRESSURE. 

The  duration  of  the  outward  flow  in  the  blast  wave  can  be 
related  to  the  peak  overpressure  (for  a  free-air  burst)  as 


(52) 


ms , 


where  tt  =  APS/1000  (in  kilopounds  per  square  inch), 
m  =  W1^  (in  cube-root  kilotons)  . 

For  a  surface  burst,  substitute  2W  for  W,  i.e. ,  m  —  1 .26m.  The  fit  is 
in  good  agreement  with  earlier  calculations  [Brode,  1959b,  1966],  but 
it  differs  significantly  from  the  1 -KT  standard  [Needham  and  Crepeau, 
1981].  The  difference  is  illustrated  in  Fig.  18  (scaled  to  1 -MT 
surface  burst).  Such  a  disparity  between  results  of  detailed 
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Figure  18.  Fit  compared  to  calculation  and  1-KT  standard  (scaled  to 
1-MT  surface  burst):  dynamic  pressure  positive  phase. 


numerical  (one-dimensional)  calculations  is  a  measure  of  the  dif¬ 
ferences  introduced  by  dissimilarities  in  boundary  and  initial  condi¬ 
tions,  equations  of  state  and  opacities,  and  by  various  treatments  of 
radiation  transport,  thermal  radiation  losses,  and  accumulated  numeri¬ 
cal  errors  in  detailed  computer  calculations. 

DYNAMIC  IMPULSE  VERSUS  PEAK  OVERPRESSURE  (FREE-AIR  BURST). 

The  dynamic  impulse  in  the  positive  phase,  defined  as 

T+D  + 

xu  H  (1/2)  f  Pu2  dt  »  (53) 

T 


can  be  approximated  as 


I 

u 


2.14AP  1 *637(m) 

_ s _ 

( 1  +  0.00434AP  1 *431 ) 

s 


psi-ms. 


(54) 


That  form  is  within  10  percent  of  the  scaled  values  from  the  detailed 
calculations  [Brode,  1959b,  1966]  for  3  <  APS  <  10,000  psi.  It  is 
high  by  nearly  20  percent  at  APS  =  100,000  psi.  Figure  19  compares 
the  dynamic  impulse  from  the  detailed  calculations  with  that  from  this 
fit  (Eq.  54). 


DYNAMIC  IMPULSE  VERSUS  SCALED  RANGE. 

A  fit  to  the  dynamic  impulse  versus  range  for  the  early  calcula¬ 
tions  [Brode,  1959b]  agrees  to  better  than  10  percent  for  0.0025  ^  r  ^ 
2  kft/KT1^.  The  relation,  when  scaled  to  a  1 -KT  free-air  burst,  is 


I  = 
u 


_ 18. 8r  _  +  92.64 

10"6  +  0.06896r3  +  0.5963r5,652  (100r)5 


2935(r  -  0 . 00597 ) (0 . 01  -  r)(0. 0003552  -  r  ) 

10~10  +  0.003377r2*5  +  155. 8r8 


m  psi-ms,  (55) 
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Figure  19.  Fit  compared  to  calculation:  dynamic  impulse  versus  peak 
overpressure  for  1-KT  free-air  burst. 


with  r  -  sr/m,  sr  in  kilofeet,  m  =  W1^,  and  W  in  kilotons.  This 
expression  is  illustrated  in  Fig.  20  and  compared  with  the  detailed 
calculation  results  to  which  it  was  fit.  The  fit  is  good  to  a  few 
percent  over  the  entire  range.  For  a  surface  burst,  m  =  (2W)1/3. 


DYNAMIC  PRESSURE  VERSUS  TIME. 

An  older  approximate  analytic  expression  for  dynamic  pressure 
versus  time  covers  the  range  2  <  APX  <  1000  psi  (0.1  <  Qs  <  3000  psi) 
[Brode,  1964]: 


Q(t)  = 


Qs(1  - 


u>)  [d  exp  (-aw)  +  (1  -  d)  exp  (-bw)]  psi, 


(56) 


in  which  co 
7  T 
't 


(t  -  T)/Dj, 

time  of  arrival  [see  Eqs.  (39)  through  (41)], 
time, 

duration  of  outward  blast  wind  [see  Eq.  (52)], 
peak  dynamic  pressure  in  pounds  per  square  inch  [see 
Eqs.  (16)  through  (18)], 


d 


a 

b 


1 


.06,0-035 


1  +  147tt‘ 


2.1  3tt 


3 


1  +  67 . 9tt 


3.5 


0.38AP  °‘8605  .  145,°*8605 
s 


K  JIAD  0.604  0.604 

5.4AP  =  3  5  Ott 

s 


f 


Equation  (56)  is  valid  for  1  Z  tt  >  0.002  ksi  (1000  2  APX  ^  2  psi). 
These  parameters  are  illustrated  in  Fig.  21  (scaled  to  1 -MT  surface 
burst ) . 

A  relatively  simple  alternative  surface  burst  formula  for  dynamic 
pressure  versus  time  can  be  derived  from  the  more  complex  fits  to 
dynamic  and  overpressure  versus  HOB,  ground  range,  and  time  [Speicher, 
1983;  Speicher  and  Brode,  1981;  Brode,  1983].  That  fit,  when 
simplified  for  zero  HOB,  becomes  essentially  that  of  the  shock  or 
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Figure  20.  Fit  compared  to  calculation:  scaled  dynamic  pressure 
impulse  versus  scaled  range  for  1-KT  free-air  burst. 
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Figure  21  Dynamic  pressure  versus  peak  overpressure  and  time 
y  (scaled  for  1-MT  surface  burst). 
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Hugoniot  form  of  the  relation  between  dynamic  pressure  and  overpres¬ 
sure,  but  extended  to  times  after  shock  arrival  [e.g.»  Eq.  (17)]: 


Q( x,  a) 


2.5AP2 

102.9  +  AP 


psi, 


(57) 


In  this  fit,  AP  is  defined  by  Eq.  (50),  except  that  the  quantity  D  in 
Eq.  (50)  should  be  replaced  by 

D'  =  D [ M .51  -  12.8x2‘5/(1  +  3.63x2*5)]  .  (58) 

This  is  not  a  wholly  accurate  form,  since  it  is  the  result  of  a 
"quick  fix"  HOB  fit  for  dynamic  pressure  time-histories,  and  has  not 
yet  been  replaced  by  a  more  exact  expression. 

MAXIMUM  TEMPERATURE  (FIREBALL  EXPOSURE)  VERSUS  PEAK  OVERPRESSURE. 

Another  parameter  of  interest  is  the  maximum  temperature  at  a 
given  range.  The  maximum  temperature  at  overpressure  levels  above 
100  psi  occurs  after  shock  arrival,  as  hotter  fireball  air  expands 
past  that  point.  That  maximum  temperature  is  somewhat  dependent  on 
yield,  since  the  fireball  of  a  megaton  burst  cools  more  slowly  than 
that  of  a  kiloton  burst  (even  if  cube-root  scaling  is  applied) .  The 
rough  fit  of  Eq.  (59)  is  appropriate  for  a  megaton  and  is  based  on 
earlier  calculations  by  Brode  [1959b]: 


0 

m 


1  090it 


3.26 


1  +  35.6-tt 


2.75 


1 03°C, 


(59) 


with  it  in  kilopounds  per  square  inch.  More  recently,  improved  opaci¬ 
ties  for  air  have  led  to  slightly  higher  fireball  temperatures  at 
corresponding  times,  so  that  temperatures  somewhat  above  those 
predicted  by  Eq.  (59)  are  likely. 

In  Fig.  22,  the  fit  is  compared  with  the  values  from  the 
radiation-hydrodynamic  calculation  [Brode,  1959b]. 
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Figure  22.  Fit  compared  to  calculation:  maximum  fireball  temperature 
versus  peak  overpressure  for  1 -MT  burst. 


The  time  to  maximum  temperature  cannot  be  rigorously  scaled;  but 
assuming  cube-root  scaling,  that  time  is  roughly  inversely  propor¬ 
tional  to  the  peak  overpressure  (free-air  burst): 

t  =  1 9 w"'  ^ / tt  =  19  m/ir  ms,  (60) 

m 

where  tt  =  APS/1000  (in  kilopounds  per  square  inch)  and  W  is  yield  in 
kilotons.  This  approximation,  which  is  based  on  a  calculation  for  a 
few  megatons,  is  less  reliable  as  the  yield  deviates  from  that. 

NEGATIVE  PHASE  UNDERPRESSURE  VERSUS  PEAK  OVERPRESSURE  AND  TIME. 

The  time-history  of  the  negative  overpressure  AP_,  or  underpres¬ 
sure,  can  be  approximated  for  a  free-air  (or  surface)  burst  by  the 
algorithm 

p  A  id  -  t) 

AP_  - - r— - -  psi,  (61) 

1  +  B  t  +  C  x 
n  n 

ambient  pressure  (14.7  psi), 

(t  -  tn)/D“, 
time  after  detonation, 
time  of  beginning  of  negative  phase, 
time  to  end  of  positive  phase, 

time  of  arrival  +  duration  of  positive  phase  =  T  +  Dp, 
duration  of  negative  phase  =  1051  m  (in 
milliseconds ) , 

wl/3  (in  cube-root  kilotons), 

0.2532APs  41 3.2( AP  /I  00) ^ 

(1  +  0.1262AP  )  +  7  “71  ~  57  ’ 

s  [l  +  668.  1  (  AP  /I  00r  1 

s 

2.481 AP 

_ s _ 

(1  +  0.004272AP  1  *7) 

s 


A  = 
n 


B  = 
n 


where  P0  = 
t  = 


t  = 


m  = 


52 


c  = 

n 


1 8. 55( AP  /100) 
s 


8 


[1  +  2.75CAP  /100) 
2 


7.335 


] 


t  =  m  1  51  .  3  + 

n  1 


2844 


AP 


0.9638 


=  T  +  D+  , 

P 

[i.e.,  Eqs.  (39)  through  (41),  plus  Eqs.  (43)  through  (43)]. 


This  fit  is  compared  with  an  early  detailed  one-dimensional  free- 
air  calculation  [Erode,  1959b]  at  overpressure  ranges  from  300  psi 
down  to  less  than  5  psi  in  Figs.  23  through  30. 


BLAST  SUMMARY. 

A  number  of  free-air  blast  parameters  (shock  radius,  time  of 
arrival,  overpressure  positive  phase,  overpressure  impulse,  dynamic 
pressure  duration,  dynamic  impulse,  and  maximum  fireball  temperature) 
are  listed  in  Table  2  for  shock  overpressures  ranging  from  0.1  to 
10,000,000  psi. 

Figure  31  summarizes  some  of  the  shock  parameters  for  a  1 -MT 
surface  burst  versus  peak  overpressure.  Several  of  the  variables 
listed  in  Tables  1  and  2  are  illustrated  in  that  figure.  They  are 
time  of  arrival  T,  in  milliseconds;  shock  radius  GR,  in  feet;  shock 
velocity  Us,  in  kilofeet  per  second;  peak  particle  velocity  us,  in 
feet  per  second;  peak  dynamic  pressure  Qs ,  in  pounds  per  square  inch; 
shock  temperature  increase  A0S,  in  degrees  centigrade;  and  maximum 
fireball  temperature  0m,  in  degrees  centigrade. 
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Table  2.  Quantities  for  1 -KT  free-air  blast  wave 


Peak 

Overpressure 

A P3  (psl) 

Shock 

Radius 
sr  ( f t )  T 

Overpressure 

Time  of  Positive 

Arrival  Phase 

(ms/KTT/3)  Dp  (ms/KT  /3) 

Overpressure 
Impulse  . /7 

I  (psi-ms  KT  J) 

Dynamic 

Pressure 

Duration 

D*  (ms/KT1/3) 

Dynamic 

Imnul3e 

Ij  (psi-ms/KT  /3) 

Maximum 

Fireball 

Temperature 

&  ( 1 0-°G) 

31 

.1 

26310 

22760 

332 

45.8 

315 

0.0274 

0.0005 

.15 

18160 

15550 

327 

56.1 

314 

0.0579 

0.0008 

.2 

14020 

11890 

322 

64.7 

313 

0.0977- 

0.0011 

.3 

9816 

8161 

315 

79.2 

311 

0.203 

0.0016 

.4 

7667 

6261 

308 

91.5 

310 

0.338 

0.002 

.6 

5467 

4317 

297 

112 

306 

0.692 

0.003 

.8 

4335 

3320 

288 

129 

302 

1.15 

0.004 

1.0 

3639 

2709 

279 

144 

298 

1.70 

0.005 

1.5 

2680 

1872 

262 

177  • 

289 

3.46 

0.007 

2 

2178 

1440 

248 

204 

280 

5.75 

0.010 

3 

1650 

994.0 

226 

249 

264 

11.8 

0.015 

4 

1370 

763.4 

209 

288 

249 

19.4 

0.019 

6 

1069 

525.9 

183 

352 

226 

38.8 

0.028 

7 

977.3 

456.3 

173 

380 

217 

49.9 

0.032 

10 

801.0 

328.9 

151 

453 

200 

87.4 

0.045 

15 

648.0 

227.1 

126 

553 

191 

156 

0.064 

20 

562.2 

174.9 

110 

637 

195 

226 

0.082 

30 

465.0 

121.5 

91.5 

778 

211 

358 

0.117 

40 

409.0 

94.05 

81.4 

895 

224 

475 

0.15 

50 

371.5 

77.22 

75.5 

998 

230 

579 

0.18 

70 

322.9 

57.50 

69.8 

1175 

231 

754 

0.24 

100 

279.8 

42.17 

67.7 

1396 

222 

961 

0.56 

150 

239.1 

29.72 

69.1 

1700 

209 

1210 

1.88 

200 

214.5 

23.21 

72.0 

1940 

202 

1390 

4.02 

300 

184.6 

16.39 

77.4 

2350 

195 

1630 

9.36 

450 

159.5 

11.57 

83.5 

2840 

192 

1790 

16.3 

700 

136.4 

7.903 

89.9 

3480 

191 

1890 

23.7 

1000 

120.4 

5.798 

94.5 

4090 

191 

1950 

29.8 

1500 

104.7 

4.064 

99.0 

4890 

191 

2060 

37.3 

2000 

94.85 

3.152 

102 

5530 

191 

2200 

43.4 

3000 

83.67 

2.106 

105 

6560 

191 

2460 

53.5 

4000 

75.06 

1.696 

107 

7370 

191 

2690 

62.0 

6000 

65.64 

1.177 

110 

8650 

192 

3060 

76.3 

8000 

59.78 

0.9091 

111 

9650 

192 

3350 

88.4 

10,000 

55.68 

0.7460 

112 

10500 

192 

3570 

99.1 

15,000 

49.17 

0.5270 

114 

12100 

192 

3970 

121 

20,000 

45.26 

0.4180 

115 

13300 

192 

4230 

141 

30,000 

40.67 

0.3100 

117 

15100 

192 

4530 

174 

40,000 

37.92 

0.2570 

118 

16400 

192 

4700 

201 

60,000 

34.12 

0.1936 

120 

18300 

192 

4880 

247 

80,000 

31.00 

0.1510 

121 

19600 

192 

4960 

286 

100,000 

28.40 

0.1206 

122 

20700 

192 

4940 

321 

150,000 

23.91 

0.07897 

124 

22500 

192 

4660 

394 

200,000 

21.22 

0.05719 

125 

23800 

192 

4300 

456 

300,000 

18.11 

0.03659 

127 

25500 

192 

3690 

561 

400,000 

16.29 

0.02583 

128 

26700 

192 

3260 

650 

600,000 

14.10 

0.01438 

129 

28200 

192 

2710 

799 

800,000 

12.76 

0.008623 

131 

29200 

192 

2390 

926 

1000,000 

11.82 

0.005415 

132 

29900 

192 

2180 

1037 

1500,000 

10.30 

0.002133 

133 

31100 

192 

1900 

1275 

2000,000 

9.352 

0.001029 

135 

31800 

192 

1760 

1477 

3000,000 

8.163 

0.0003522 

136 

32750 

192 

1670 

1816 

4000,000 

7.414 

0.0001623 

138 

33300 

192 

1660 

2103 

6000,000 

6.475 

0.0000514 

140 

34100 

192 

1750 

2587 

8000,000 

5.882 

0.00002485 

141 

34500 

192 

1920 

2996 

10,000,000 

5.460 

0.00001360 

142 

34800 

192 

2140 

3357 

Mote:  Shock  radius,  see  Eq.  (36);  time  of  arrival,  see  Eq.  ( 4 1 ) ;  overpressure  positive  phase, 

see  Eq.  (43);  overpressure  impulse,  see  Eq.  (48);  dynamic  pressure  duration,  see  Eq.  (52); 
dynamic  impulse,  see  Eq.  (55);  maximum  fireball  temperature,  see  Eqs.  (20)  and  (24),  or  Eq.  (59). 
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Fi qure  31  .  Blcist  psrdrriGtGrs  for  1  MT  surfdCG  burst. 
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SECTION  4 

BLAST  PARAMETERS  AS  FUNCTIONS  OF  BURST  HEIGHT 


PEAK  OVERPRESSURE  VERSUS  SCALED  BURST  HEIGHT  AND  SCALED  RANGE. 

An  earlier  fit  to  the  peak  overpressure  HOB  contours  was  valid  up 
to  pressures  of  10,000  psi  [Speicher  and  Brode,  1980b].  It  was  sub¬ 
sequently  extended  to  much  higher  pressures  [Speicher  and  Brode, 
1984a].  The  latter  fit  was  illustrated  to  3,000,000  psi,  but  the 
basis  for  that  fit  was  restricted,  at  the  time,  to  the  old  data  (valid 
to  10,000  psi),  to  a  calculation  near  zero  burst  height,  which  is 
called  BM-3  [Pyatt,  1983],  and  to  some  low-altitude  HULL  calculations 
at  25  and  50  ft  scaled  height  of  burst  (SHOB)  [Pyatt  and  Wilkins, 
1983].  Interpolation  for  bursts  between  zero  and  200-ft  SHOB  was 
modestly  successful,  as  borne  out  by  comparison  with  results  of  sub¬ 
sequent  calculations  [Fry,  Kamath,  and  Book  1985;  Fry,  1986].  Those 
results,  however,  were  not  in  as  good  agreement  with  some  later  HULL 
calculations  at  the  low-burst  altitudes  [Pyatt,  1985].  Figures  32 
through  40  span  the  range  of  the  current  fit,  showing  contours  of 
constant  peak  overpressure  plotted  against  scaled  burst  height  and 
scaled  ground  range.  Comparisons  with  calculation  results  from  the 
Naval  Research  Laboratory  (NRL) ,  Washington,  D.C.,  and  S-Cubed,  La 
Jolla,  California,  are  plotted  in  Figs.  32  through  36.  In  those 
figures,  the  lines  that  indicate  the  onset  of  Mach  reflection  (Xm)  and 
the  locus  of  points  where  the  second  peak  equals  the  first  peak  (Xe) 
are  sometimes  shown.  Those  curves  are  approximated  analytically 
(along  with  the  time  dependence  of  pressure)  in  Eq.  (63).  For  NRL 
calculations,  agreement,  in  most  cases,  is  within  20  percent  on  peak 
overpressure. 

The  fit  for  peak  overpressure  takes  the  form: 

psi, 

(62) 


AP 


10.47 

a(z) 


b(z) 

c(z) 


d(z)  x  e(z) 


1  +  f(z)  x  r 


g(z) 


h(z, 


y) 


J 

nk(y) 
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□  NRL  40,000  psi 
0  NRL  106  psi 


d  70,000  psi 
O  S3  42,000  psi 
(?S3  106  psi 


^002  0.03  0.04 

Scaled  ground  range  (kft/KT  '  ) 

Figure  32.  Ultra-high  peak  overpressure  HOB  versus  ground  range  contours, 
scaled  to  1  KT. 
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□  NRL  40,000  psi 
0  S3  25,000  psi 
O  S3  42,000  psi 


Figure  33.  Extremely  high  peak  overpressure  HOB  versus  ground  range  contours 
scaled  to  1  KT. 
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Figure  36.  Intermediate  high  peak  overpressure  burst  height  versus  ground 
range  contours,  scaled  to  1  KT. 
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Figure  37.  Intermediate  peak  overpressure  burst  height  versus  ground  range 
contours,  scaled  to  1  KT. 
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Figure  34.  Very  high  peak  overpressure  HOB  versus  ground  range  contours, 
scaled  to  1  KT. 
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Figure  35.  High  peak  overpressure  HOB  versus  ground  range  contours, 
scaled  to  1  KT. 
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Figure  38.  Intermediate  low  peak  overpressure  burst  height  versus  ground 
range  contours,  scaled  to  1  KT. 
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Figure  39.  Low  peak  overpressure  burst  height  versus  ground  range  contours, 
scaled  to  1  KT. 
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Figure  40.  Very  low  peak  overpressure  burst  height  versus  ground  range  contours, 
scaled  to  1  KT. 
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where  r 


s  scaled  slant  range  in  kilofeet  per  cube-root  kiloton 
=  /x^  +  y2  f 

x  =  scaled  ground  range  in  kilofeet  per  cube-root  kiloton,  or 
GR/m/1 000 , 

y  =  scaled  burst  height  in  kilofeet  per  cube-root  kiloton,  or 
H/m/1000, 

m  =  W1^  in  cube-root  kilotons  (the  scale  factor), 

W  s  yield  in  kilotons, 

GR  =  ground  range  in  feet, 

H  ^  burst  height  in  feet, 
z  3  H/GR  =  y/x, 

and  where 


a(z)  =  1.22- 
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h(z,  r,  y)  = 


8 . 808z  _  0.2905  +  64.67z^  1 . 3  8  9  z 


1  +  154. 5z 
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OVERPRESSURE  VERSUS  TIME,  SCALED  RANGE,  AND  BURST  HEIGHT. 

An  earlier  time-history  fit,  published  by  Speicher  and  Brode 
[1981],  included  the  double  peaks  in  the  high  overpressure  reflection 
region  of  the  HOB-range  map  (Fig.  41),  but  was  limited  to  peak  over¬ 
pressures  less  than  10,000  psi  and  more  than  1  psi  (or  to  0.1  psi  for 
a  surface  burst).  Four  sources  formed  the  basis  of  that  fit:  (1) 
results  of  a  calculation  at  200  ft  SH0B  [McNamara,  Jordano,  and  Lewis, 
1977],  (2)  a  composite  free-air  blast,  1 -KT  "standard"  [Needham, 
Havens,  and  Knauth,  1975],  (3)  the  results  of  a  large-charge,  high- 
explosive  (HE)  test  series  called  DIPOLE  WEST  [Reisler,  1980],  and  (4) 
tests  with  8-lb  spheres  of  PBX-9404  (plastic  bonded  explosive)  con¬ 
ducted  and  analyzed  by  Carpenter  [1976],  which  concentrated  on  the 
higher  overpressure  regions  of  the  early  Mach  reflections  (100  to 
2500  psi ) . 

The  fit  reported  here  includes  a  recent  extension  to  much  higher 
overpressures  [Speicher  and  Brode,  1984b].  It  is  unlikely  that  the 
fit  will  be  useful  at  overpressures  higher  than  200,000  psi,  because 
radiation  transport  is  a  dominant  factor  that  near  the  burst.  In  a 
1-MT  surface  burst,  radiation  dominates  the  fireball  expansion  out  to 
a  radius  of  about  250  ft  (until  the  pressure  drops  to  approximately 
200,000  psi).  The  shock  wave  is  not  fully  developed  prior  to  that. 

*In  Fig.  41 ,  the  line  depicting  the  locus  of  points  where  the  sec¬ 
ond  peak  is  just  equal  to  the  first  peak  (Xe)  was  first  drawn  by 
Carpenter  [1976] . 
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Figure  41.  High  peak  overpressures  versus  scaled  burst  height  and  scaled  ground 
range  (ideal  surface). 
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Although  reflected  pressures  can  rise  above  the  incident  shock  pres¬ 
sure,  much  of  the  energy  of  an  incident  shock  of  200,000  psi,  on 
reflection,  radiates  away,  allowing  the  density  to  rise  without  sig¬ 
nificantly  altering  reflected  overpressure. 

The  changes  to  the  fit  given  here  were  modeled  after  the  close-in 
results  of  the  surface  burst  BM~3  calculation  [Pyatt,  1983]  and  HULL 
calculations  at  25  and  50  ft  SHOB  [Pyatt  and  Wilkins,  1983].  However, 
the  fit  does  not  correspond  well  with  calculations  of  the  early  Mach 
reflection;  it  predicts  too  great  an  impulse  around  the  second  peak. 

At  an  SHOB  of  25  ft  and  scaled  ground  ranges  of  30  to  100  ft,  the 
partial  impulse  to  2  ms  at  1  KT  is  high  by  as  much  as  a  factor  of  2. 
That  partial  impulse  correlates  closely  with  blast  damage  to  hardened 
surface-flush  structures  such  as  missile  silos,  so  the  error  can  have 
serious  consequences.  This  section  of  the  fit  is  being  modified,  and 
copies  of  the  improved  form  will  be  distributed  later.  Elsewhere  (at 
the  surface  or  at  SHOBs  above  100  ft  and  scaled  ground  ranges  beyond 
100  ft),  the  fit  is  still  quite  good. 

The  time-history  of  overpressure  at  various  ranges  (X)  and  burst 
heights  (Y)  can  be  expressed  as 


AP(X,  Y,  t) 


1/3 

AP  (1  +  a) (bv  +  c)  for  X  >  Xm  and  Y  <  380  ft/KT  J  , 
s 

1/3 

aP  x  b  for  X  <  Xm  or  Y  >  380  ft/KT 

(63) 
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X 

GR 

m 
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overpressure  in  pounds  per  square  inch, 

scaled  ground  range  in  feet  per  cube-root  kiloton, 

or  GR/m, 

ground  range  in  feet, 

W1/3  in  cube-root  kilotons, 
yield  in  kilotons, 
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Y  =  scaled  burst  height  in  feet  per  cube-root  kiloton,  or 
H/m, 

H  =  height  of  burst  in  feet, 
t  =  time  after  detonation  in  milliseconds, 

APS  =  peak  overpressure  in  pounds  per  square  inch 
[see  Eq.  (62)], 
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scaled,  in  feet  per  cube-root  kiloton, 
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t  =  scaled  time  of  arrival  in  milliseconds  per  cube-root 


kiloton,  based  on  Eq.  (41), 

=  u(r)  for  X  <  Xm, 
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In  the  above  definitions  for  Eq.  (63),  the  expression  for  free-air 
time  of  arrival  u  is  the  same  as  Eq.  (42),  and  the  scaling  factor  for 
it  is  m  =  1 ;  the  expression  w  is  also  the  same  as  Eq.  (42),  but  with 
m  =  2^3. 

Examples  showing  how  successful  the  fit  is  for  Eqs.  (62)  and  (63) 
are  provided  in  Figs.  42  through  87.  In  those  figures,  the  fit  is 
compared  with  six  sources  of  "data":  the  S-Cubed  [Pyatt,  1983] 
surface-burst  calculation  (BM-3)  scaled  to  1  KT  (Figs.  42  through  49); 
the  DNA  1 -KT  standard  [Needham  and  Crepeau,  1981]  (Figs.  50  through 
58);  the  25  and  50  ft  SHOB  S-Cubed  [Pyatt,  1983]  calculations 
(Figs.  59  through  70);  the  tests  of  PBX-9404  8-lb  HE  spheres  [Car¬ 
penter,  1976]  (Figs.  71  and  72);  the  HULL  calculation  for  60-m  SHOB 
[McNamara,  Jordano,  and  Lewis,  1977]  (Figs.  73  and  74);  and  the  DIPOLE 
WEST  data  [Reisler,  1980]  (Figs.  75  and  76). 

The  plots  display  overpressure  versus  time  and  impulse  versus 
time.  The  diversity  of  data  sources  and  the  extensive  range  of  the 
fit  mitigate  against  its  being  everywhere  close.  However,  the  in¬ 
herent  uncertainty  and  expected  variation  in  actual  pressure  histories 
and  impulses  are  far  larger  than  the  10  to  20  percent  disparities 
between  this  analytic  fit  and  the  calculation  results  shown  in  these 
comparison  plots.  (Average  discrepancies  between  fit  and  calculated 
pressures  are  of  the  order  of  only  4  percent). 

The  simple  fit  for  positive  phase  overpressure  impulse  in  a 
surface  burst  [Eq.  (48)]  predicts  somewhat  higher  values  than  those 
for  the  integral  of  the  time-history  fit  for  zero  burst  height  [time 
integral  of  Eqs.  (62)  and  (63)].  The  latter,  while  much  more  complex, 
was  fit  to  the  BM-3  cratering  calculation  impulses  [Pyatt,  1983]  in 
the  close-in  range  (scaled  ground  range  between  25  and  100  ft,  or  peak 
pressure  between  1500  and  150,000  psi)  .  The  lower  peak  overpressure 
in  the  detailed  calculations  may  not  be  correct,  although  further 
calculations  are  being  made  to  verify  that  region  of  the  BM~3 
calculation.  It  has  been  suggested  that  the  lower  impulses  may  be  due 
to  the  quenching  effect  of  cratered  material  ejected  into  the 
fireball. 
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Overpressure  (psi) 


Figure  42.  Fit  compared  to  S-Cubed  calculation  (scaled):  overpressure  and 
scaled  partial  impulse  versus  scaled  time  to  0.07  ms,  for  aP  « 
548,300  psi,  y  =  0,  x  =  18.896  ft/KT1'3. 
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Overpressure  (psi) 


Figure  43.  Fit  compared  to  S-Cubed  calculation  (scaled):  overpressure  and 
scaled  partial  impulse  versus  scaled  time  to  3.5  ms,  for  aP  ~ 
548,300  psi,  y  =  0,  x  =  18.896  ft/KT1/3. 
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Scaled  impulse  (psi-ms/KT 


Overpressure  (psi)  and  scaled  impulse  (psi-ms/KT  ) 


Figure  44. 


Fit  for  overpressure  and  scaled  partial 
to  70  ms,  for  APS  ~  548,300  psi,  y  =  0, 


impulse  versus  scaled  time 
x  =  18.896  ft/KT1/3. 
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2  3  4  5 

Scaled  time  (ms/KT1^) 


Figure  46.  Fit  compared  to  S-Cubed  calculation  (scaled):  overpressure  and 
scaled  impulse  versus  scaled  time  to  7  ms,  for  ~  99,740  psi, 
y  =  0,  x  =  33.068  ft/KT1/3. 
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Overpressure  (psi)  and  scaled  impulse  (psi-ms 
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Figure  47.  Fit  for  overpressure  and  scaled  impulse  versus  scaled  time  to 
140  ms,  for  aPs  ^  99,740  psi,  y  =  0,  x  =  33.068  ft/KT1/'3. 
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Overpressure  (psi)  and  scaled  impulse  (psi-ms/KT1' 


2  3  4  5 

Scaled  time  (ms/KT1^) 


Figure  48. 


Fit  compared  to  S-Cubed  calculation: 
versus  scaled  time  to  7  ms,  for  AP$  ~ 


overpressure  and  scaled  impulse 
10,760  psi,  y  =  0,  x  =  69  ft/KT^/^. 
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Overpressure  (psi)  and  scaled  impulse  (psi-ms/KT  /J) 


Scaled  time  (ms/KT^) 

Figure  49.  Fit  for  overpressure  and  scaled  impulse  versus  scaled  time  to  70  ms, 
for  aPs  ~  10,760  psi,  y  =  0,  x  =  69  ft/KT1/3. 
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SHOB  =  0  kft,  SGR  =  0.362  kft 


Figure  51. 


Fit  compared  to  DNA  1-KT  standard  (2W):  overpressure  versus  scaled 
time  to  100  ms,  for  AP$  ^  98.59  psi. 
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Figure  52.  Fit  compared  to  DNA  1-KT  standard  (2W):  scaled  impulse  versus  scaled 
time  to  100  ms,  for  aP$  ~  98 . 59  psi. 


SHOB  =  0  kft,  SGR  =  0.362  kft 
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SHOB  =  0  kft,  SGR  =  1.03  kft 


Figure  53. 


Fit  compared  to  DNA  1-KT  standard  (2W):  overpressure  versus  scaled 
time  to  250  ms,  for  ~  10.01  psi. 
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SHOB  -  0  kft,  SGR  =  1.03  kft 


Figure  54. 


Fit  compared  to  DNA  1-KT  standard  (2W):  scaled  impulse  versus 
scaled  time  to  250  ms,  for  aP$  ~  10.01  psi . 
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Overpressure  (psi) 


SHOB  =  0  kft,  SGR  =  4.554  kft 


Figure  55.  Fit  compared  to  DNA  1-KT  standard  (2W):  overpressure  versus 
scaled  time  to  450  ms,  for  AP$  ~  1.036  psi. 
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Scaled  impulse  (psi-ms) 


240 


SHOB  =  0  kft,  SGR  =  4.554  kft 


Figure  56. 


Fit  compared  to  DNA  1-KT  standard  (2W):  scaled  impulse  versus 
scaled  time  to  450  ms,  for  aP^  ~  1.036  psi . 
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Overpressure  (psi) 


SHOB  =  0  kft,  SGR  =  27  kft 


Figure  57. 


Fit  compared  to  DNA  1-KT  standard  (2W):  overpressure  versus 
scaled  time  to  500  ms,  for  AP$  0.1086  psi. 
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Scaled  impulse  (psi-ms) 


30 


SHOB  =  0  kft,  SGR  =  27  kft 


Figure  58. 


Fit  compared  to  DNA  1-KT  standard  (2W):  scaled  impulse  versus 
scaled  time  to  500  ms,  for  aP$  ~  0.1086  psi . 
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Overpressure  (ksi) 


Overpressure  and  impulse  vs  time  (SHOB  =  25  ft,  SGR  =  6.56  ft) 


Scaled  time  (ms) 


Figure  59.  Fit  compared  to  scaled  SAI  calculation:  overpressure  and 
scaled  impulse  versus  scaled  time  to  .13  ms,  for  25  ft  HOB 
AP  %  1 ,785,000  psi. 
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Scaled  impulse  (ksi-ms) 


Overpressure  and  impulse  vs  time  (SHOB  =  25  ft,  SGR  =  6.56  ft) 


Figure  60.  Fit  compared  to  SAI  calculation:  overpressure  and 

scaled  impulse  versus  scaled  time  to  .41  ms,  for  25  Ft 

SHOB,  AP  %  1,785,000  psi. 
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Scaled  impulse  (ksi-ms) 


Overpressure  (psi) 


Overpressure  and  impulse  vs  time  (SHOB  =  25  ft,  SGR  -  32.82  ft) 
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Figure  61 . 


Fit  compared  to  SAI  calculation:  overpressure  and 

scaled  impulse  versus  scaled  time  to  .55  ms,  for  25  ft 

SHOB,  AP  jfc  233,300  psi. 
s 
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Scaled  impulse  (psi-ms) 


Overpressure  (psi) 


Overpressure  and  impulse  vs  time  (SHOB  =  25  ft,  SGR  =  32.82  ft) 


Figure  62.  Fit  compared  to  SAI  calculation:  overpressure  and 

scaled  impulse  versus  scaled  time  to  1.6  ms,  for  25  Ft 
SHOB,  APs  £  233,300  psi. 
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Scaled  impulse  (psi-ms) 


Overpressure  (ksi) 


Overpressure  vs  time-TOA  (SHOB  =  50  ft,  SGR  =  82.02  ft) 


Scaled  time-TOA  (ms) 


Figure  63.  Fit  compared  to  S-CUBED  calculation:  overpressure  versus 
scaled  time  minus  time-of-arrival  to  7  ms,  for  50  ft  SHOB, 
APg  ^  9,000  psi. 
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Scaled  impulse  (ksi-ms) 


Scaled  impulse  vs  time-TOA  (SHOB  =  50  ft,  SGR  =  82.02  ft) 


Figure  64.  Fit  compared  to  S— CUBED  calculation:  scaled  impulse 
versus  scaled  time  minus  time-of-arrival  to  7  ms ,  for 
50  ft  SHOB,  AP  %  9,000  psi. 

o 
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Scaled  impulse  (ksi-ms) 


Scaled  impulse  vs  time-TOA  (SHOB  =  50  ft,  SGR  -  82.02  ft) 


Figure  65.  Fit  compared  to  S-CUBED  calculation:  scaled  impulse 

versus  scaled  time  minus  time-of-arrival  to  100  ms, 

for  50  ft  SHOB,  AP  %  9,000  psi. 
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Overpressure  (psi) 


SHOB  =  0.107  kft,  SGR  =  0.151  kft 


Figure  66.  Fit  compared  to  8-lb  charge  data:  overpressure  versus  scaled  time 
to  10  ms,  for  107-ft  SHOB,  aP$  ~  1199  psi. 
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Scaled  impulse  (psi-ms) 


6000 


SHOB  =  0.107  kft,  SGR  =  0.151  kft 


Figure  67.  Fit  compared  to  8-lb  charge  data:  scaled  impulse  versus  scaled  time 
to  10  ms,  for  107-ft  SHOB,  aP$  «  1199  psi . 
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Overpressure  (psi) 


Figure 


58 


Fit  compared  to  General  Electric-TEMPO  calculation:  overpressure 


versus 


scaled  time  to  70  ms,  for  60-m  SHOB,  aP  ~  110.5  psi 
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Scaled  impulse  (psi-ms) 


2400 


SHOB  =  0.1969  kft,  SGR  =  0.3281  kft 


Scaled  time  (ms/KT^) 


Figure  69*  Fit  compared  to  General  Electric-TEMPO  calculation:  scaled  impulse 
versus  scaled  time  to  70  ms,  for  60-m  SHOB,  aPs  %  110.5  psi. 
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SHOB  =  0.15  kft,  SGR  =  1.062  kft 


Figure  70 . 


Fit  compared  to  DIPOLE  WEST  data:  overpressure  versus  scaled  time 
to  220  ms,  for  150-ft  SHOB,  AP$  ~  10.15  psi. 
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Scaled  impulse  (psi-ms) 


900 


SHOB  =  0.15  kft,  SGR  =  1.062  kft 


240 


Figure  71  .  Fit  compared  to  DIPOLE  WEST  data:  scaled  impulse  versus  scaled 
time  to  220  ms,  for  150-ft  SHOB,  aP$  %  10.15  psi. 
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Overpressure  (psi) 


Overpressure  (200-ft  SHOB) 


Figure  72.  Fit  compared  to  200-ft  SHOB  Kaman  AviDyne  calculation:  peak 
overpressure  (from  1600  psi)  versus  scaled  ground  range 
(0  to  0.3  kft) . 
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Overpressure  (200-ft  SHOB) 


Figure  73. 


Fit  compared  to  200-ft  SHOB  Kaman  AviDyne  calculation:  peak 
overpressure  (100  to  15  psi)  versus  scaled  ground  range  (0.3  to  0.9  kft) 
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Figure 


Overpressure  (200-ft  SHOB) 


74  Fit  compared  to  200-ft  SHOB  Kaman  AviDyne  calculation:  peak 
overpressure  (10  to  1.5  psi)  versus  scaled  ground  range 

(1  to  3.7  kft) . 


Ill 


Scaled  overpressure  impulse  (200-ft  SHOB) 


Figure  75  .  Fit  compared  to  200-ft  SHOB  Kaman  AviDyne  calculation:  scaled 
overpressure  impulse  versus  scaled  ground  range  (to  1  kft). 


Scaled  overpressure  impulse  (psi-s/KT 


Scaled  overpressure  impulse  (200-ft  SHOB) 


Figure  76.  Fit  compared  to  200-ft  SHOB  KamanAvi Dyne  calculation:  impulse 
versus  scaled  ground  range  (1  to  i.b  kit). 
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e/l±>i/ij>i)  goHS 


Impulse  comparison  at  T  =  0.5  ms 


Scaled  ground  range  (kft/KT1^) 

- Analytic  [Eqs.  (62)  and  (63)] 

+■  ■+*  Carpenter  [1976] 


Figure  77.  HOB  .i soimpul se  curves:  fit  compared  to  8-lb  charge  data,  showing 
partial  impulse  contours  for  0.5  ms/KTl/3. 
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SHOB  (kft/KT 


Impulse  comparison  at  T  =  2  ms 


Scaled  ground  range  (kft/KT1^) 

- Analytic  [Eqs.  (62)  and  (63)] 

4- 4- Carpenter  [1976] 


Figure  78.  HOB  isoimpulse  curves:  fit  compared  to  8-lb  charge  data,  showing 
partial  impulse  contours  for  2  ms/KT1/^. 
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- - Analytic  [Eqs.  (62)  and  (63)] 

-h  +  Carpenter  [1976] 


Figure  79.  HOB  isoimpulse  curves:  fit  compared  to  8-lb  charge  data,  showing 
partial  impulse  contours  for  10  ms/KT^/3. 
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Partial  impulse  at  T  =  0.5  ms 


Scaled  ground  range  (kft/KT1^) 


Figure  80  .  Partial  impulse  contours  at  very  high  overpressure  for  0.5  ms/KT^-/^. 
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SHOB  (kft/KT ,/J) 


Partial  impulse  at  T  =  2  ms 


Scaled  ground  range  (kft/KT1^) 


Figure-  81  . 


Partial  impulse  contours  at  very  high  overpressure  for  2  ms/KT*/3. 
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Fi gure  82 . 


Partial  impulse  contours  at  very  high  overpressure  for  10  ms/KT^/8. 
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Although  the  Kaman  AviDyne  (KA)  calculations  by  Smiley,  Tomayko, 
and  Ruetenik  [1982a-d]  and  the  NRL  calculations  of  Fry,  Kamath,  and 
Book  [1985]  were  not  available  when  this  fit  was  generated,  com¬ 
parisons  of  the  fit  to  those  calculations  are  quite  favorable. 

Figures  77  through  81  show  the  peak  overpressures  and  total  positive 
overpressure  impulse  as  functions  of  scaled  ground  range  for  the  KA 
calculations  of  the  fit  [Eqs.  (62)  and  (63)]  at  an  SHOB  of 
200  ft/KTly/3.  The  largest  disparity  is  shown  in  Fig.  77  at  a  range 
just  after  the  onset  of  Mach  reflection,  where  the  fit  (based  on  the 
8-lb  charge  data)  shows  considerable  pressure  increase,  while  the  KA 
calculations  for  the  same  range  do  not.  The  KA  calculations  may 
suffer  most  from  lack  of  fine  grid  zoning  in  this  region. 

The  overpressure  impulse  delivered  with  a  time  interval  com¬ 
parable  to  the  time  to  fail  for  a  structure  is  relevant  to  dynamic 
response  analyses.  For  stiff  structures  designed  to  survive  at  high 
overpressures,  the  scaled  times  of  interest  are  the  first  few 
milliseconds.  Figures  82  through  87  illustrate  scaled  partial  impulse 
contours  (versus  SHOB  and  scaled  ground  range)  for  0.5,  2,  and 
10  ms/KT1^.  Figures  82  through  84  compare  the  integral  of  Eqs.  (62) 
and  (63)  with  the  partial  impulses  scaled  from  the  8-lb  charge  tests 
using  PBX-9404  [Carpenter,  1976].  Figures  85  through  87  present 
partial  impulse  contours  to  higher  overpressure  levels  (at  closer 
ranges ) . 

Appendix  A  contains  test  values  of  this  fit  for  the  purpose  of 
code  checking.  A  Fortran  subroutine  for  the  same  fit  is  given  in 
Appendix  B. 

PEAK  DYNAMIC  PRESSURE  VERSUS  SCALED  BURST  HEIGHT  AND  SCALED  RANGE. 

Speicher  [1983]  has  provided  an  improved  fit  to  the  horizontal 
peak  dynamic  pressure  from  calculations  by  KA  [Smiley,  Tomayko,  and 
Ruetenik,  1982;  Smiley,  Ruetenik,  and  Tomayko,  1984a, b];  and  by 
[Brode,  1959b].  It  follows  the  general  form  suggested  earlier  by 
Brode  [1983].  However,  it  extends  the  fit  over  the  entire  range 
inside  the  regular  reflection  region,  and  is  valid  for  dynamic  pres¬ 
sures  from  1000  psi  to  0.05  psi.  The  form  for  the  peak  dynamic  pres¬ 
sure  Qg  is  as  follows: 
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Ar^  C 

Q  =  Q  =  - —  +  —  f  n  r  x  >  X 

3  W1  .  „  E  F  1  ■  V 

1  +  Br  r  H 


<3  =  Q  exp 

s  m  K 


GL' 


1  +  649L1  1  +  HLJ 


4.01LJ  ^  _  -6 

+  7.67x10 


K  +  L 


3.22  K 


for  x  <  X  , 

q 


(64) 


where  r 


x 

y 


scaled  slant  range  in  kilofeet  per  cube-root  kiloton, 

yx2  +  y2 ' 

scaled  ground  range  in  kilofeet  per  cube-root  kiloton, 
scaled  burst  height  in  kilofeet  per  cube-root  kiloton, 
approximate  interface  between  regular  and  Mach  reflection 
in  kilofeet  per  cube-root  kiloton, 


63. 5y 


7.26 


1  +  67.1 1 y 


4.746 


+  0.6953y 


0.808 


Qm  =  Q-|  evaluated  at  x  =  Xq  (at  M  =  1), 

i  . e.  ,  for  r  =  /Xq  +  y2 , 

M  =  Xq/x, 

L  =  logio(M), 


A  =  -236.1  + 


17.72M 


0.593 


1  +  10. 4M 


3.124  ’ 


B  =  12.27 


21 . 69M 


2.24 


1  +  6.976M 


0.484 


9 


C  =  20.26  + 


14. 7M 


1  +  0.08747M 


3.05  ’ 
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D  =  -1.137 


0.5606M0,895 

1  +  3.046M7,48 


E  =  1  .731  + 


1 0.84M 


1.12 


1  +  1 2.26M 


0.0014  ’ 


F 


2.84  + 


0.855M 


0.9 


1  +  1 . 05M 


2.84 


G  =  50  - 


I843y 


2.153 


1  +  3.95y 


5.08  ’ 


H  =  0.294  +  - -----E 


8.7 


1  +  115. 3y 


6.186  ’ 


I  = 


-3.32  4  + 


4.77 


1+211  .8y 


5. 1  66 


J  =  1  .955  + 


169. 7y 


9.317 


1  +  97 . 36y 


6.513  ’ 


K  =  8.123x10 


-6  0.001 61 3y 


6.428 


1  +  60.26y 


7.358  * 


In  the  above  equation,  Xq  is  within  16  percent  of  Xm  [as  Xm  is  defined 
in  Eq.  (63)]. 

The  closeness  of  this  fit  is  illustrated  in  Figs.  88  through  106. 
In  those  figures,  the  data  points  are  represented  by  triangles  (A). 
Figures  88  through  90  compare  the  peak  dynamic  pressure  versus  range 
for  a  surface  burst  with  the  detailed  calculation  results  of  KA 
[Ruetenik,  1984]  and  the  early  calculation  by  Brode  [1959b],  The  fit 
is  very  close  at  all  ranges  from  190  to  3000  ft/KT^^  i.e.,  from  1400 
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PEAK  DYNAMIC  PR.  -  RANGE  C SHC8*a  FT! 


Source:  Erode  [1959b]  and  Ruetenik  [1984]. 

Figure  83.  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  surface  burst,  high-pressure  region. 


PEAK  DYNAMIC  PR.  -  RANGE  ( SHOB-3  FT) 


Note:  a  =  calculation  result. 

Source:  Brode  [1959b]  and  Ruetenik  [1984]. 

Figure  84-  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation 
for  surface  burst,  intermediate  region. 
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PEAK  OYNAHIC  PR.  -  RANGE  C  SHO6=0  FT) 


SCALED  GR  (FT) 

Note:  a  =  calculation  result. 

Source:  Brode  [1959b]  and  Ruetenik  [1984]. 


Figure  85.  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  surface  burst,  low-pressure  region. 


PEAK  QYNADIC  PR.  -  RANGE  ( SHO6-200  FT) 


Note:  a  =  calculation  result. 

Source:  Smiley,  Tomayko,  and  Ruetenik  [1982  a-d] . 


Figure  86  .  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  scaled  burst  height  of  200  ft/KT1/3. 
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Note:  a  =  calculation  result. 

Source:  Smiley,  Tomayko,  and  Ruetenik  [1982  a-d] . 

Figure  37.  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  scaled  burst  height  of  250  ft/KT1/3. 


PEAK  OTNAfllC  PR.  -  RANGE  C  SHO6-300  m 


Note:  a  =  calculation  result. 

Source:  Smiley,  Tomayko,  and  Ruetenik  [1982  a*d] . 


Figure  88.  Peak  dynamic  pressure  fit  compared  to  Kaman  Avidyne  calculation, 
for  scaled  burst  height  of  300  ft/KT^-/3. 
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PEAK  DYNAMIC  PR.  -  RANGE  C  SHO&n400  FT) 


Note:  a  =  calculation  result. 

Source:  Smiley,  Tomayko,  and  Ruetenik  [1982a-d] . 


Figure  89.  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  scaled  burst  height  of  400  ft/KT^/3. 


PEAK  DYNAMIC  PR.  -  RANGE  (  SHOB-600  FT) 


Note:  a  =  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [1984b] . 


Figure  90.  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  scaled  burst  height  of  500  ft/KT1/3. 
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SCALED  OR  (FT) 

Note:  a  =  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [  1984b] . 

Figure  91.  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  scaled  burst  height  of  600  ft/KT1/'3. 


PEAK  OTNAflIC  PR.  -  RANGE  (  SHO&-700  FT) 


Note:  a  =  calculation  result. 

Source:  Smiley,  Tomayko,  and  Ruetenik  [1982a-d] . 


Fiqure  92.  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  scaled  burst  height  of  700  ft/KT  /  . 
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PEAK  0TNAT1IC  Pf?.  -  RANGE  (  SHO&=750  FT) 


Note:  a  =  calculation  result. 


Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a] . 


Figure  93, .  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  scaled  burst  height  of  750  ft/KT^. 


PEAK  OYNAfllC  PR.  -  RANGE  CSHO&-1000  FT) 


SCALED  GR  CFT1 

Note:  a=  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a] . 


Figure  94.  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  scaled  burst  height  of  1000  ft/KT^. 
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Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a] . 

Figure  95-  Peak  dynamic  pressure  fit  compared  to  Kaman  AviDyne  calculation, 
for  scaled  burst  height  of  1250  ft/Kfl/3. 


CONTOURS-  1000.  500.  200.  100  PS I 


SCALED  GR  (  cm 

Note:  a  =*  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a,  b] ;  Smiley,  Tomayko, 
and  Ruetenik  [1982  a-d] ;  Ruetenik  [1984];  and  Brode  [1959b]. 


Figure  9 6.  HOB  isopicnic  contours  from  Kaman  AviDyne  calculations  compared 
to  fit  for  peak  dynamic  pressure  (100  to  1000  psi). 
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CONTOURS-  100.  50.  30.  20  PS I 


SCALED  OR  (  KFT3 


Note:  a  =  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a,  b] ;  Smiley,  Tomayko, 
and  Ruetenik  [1982  a-d] ;  Ruetenik  [1984];  and  Brode  [1959b] . 


Figure  97.  HOB  isopicnic  contours  from  Kaman  AviDyne  calculations  compared 
to  fit  for  peak  dynamic  pressure  (20  to  100  psi). 


CONTOURS-  20.  16.  10.  8  PSI 


SCALED  OR  (  Kim 

Note:  a  =  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a,  b] ;  Smiley,  Tomayko, 
and  Ruetenik  [1982  a-d] ;  Ruetenik  [1984] ;  and  Brode  [1959b] . 

Figure  98.  HOB  isopicnic  contours  from  Kaman  AviDyne  calculations  compared 
to  fit  for  peak  dynamic  pressure  (8  to  20  psi). 
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Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a,  b] ;  Smiley,  Tomayko, 
and  Ruetenik  [1982  a-d] ;  Ruetenik  [1984] ;  and  Brode  [1959b] . 

Figure  99.  HOB  isopicnic  contours  from  Kaman  AviDyne  calculations  compared 
to  fit  for  peak  dynamic  pressure  (3  to  8  psi). 


CONTOURS-  3.  2.  1.5,  1  PSI 


Note:  a  =  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a,  b] ;  Smiley,  Tomayko, 
and  Ruetenik  [  1982  a-d] ;  Ruetenik  [1984];  and  Brode  [1959b]. 


Figure  100.  HOB  isopicnic  contours  from  Kaman  AviDyne  calculations  compared 
to  fit  for  peak  dynamic  pressure  (1  to  3  psi). 
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CONTOURS-  1.  .7.  .5.  .4  PSI 


SCALED  GR  (  KFD 

Note:  a  =  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a,  b] ;  Smiley,  Tomayko, 
and  Ruetenik  [1982  a-d] ;  Ruetenik  [1984] ;  and  Brode  [1959b] . 


Figure  101.  HOB  isopicnic  contours  from  Kaman  AviDyne  calculations  compared 
to  fit  for  peak  dynamic  pressure  (0.4  to  1  psi). 
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down  to  0.05  psi.  Figures  91  through  100  compare  the  peak  dynamic 
pressure  fit  for  various  scaled  burst  heights  from  200  to 
1250  ft/KT1/3,  corresponding  to  KA  calculations  [Smiley,  Tomayko,  and 
Ruetenik,  1982;  Smiley,  Ruetenik,  and  Tomayko,  I984afb].  Again,  the 
fit  is  very  good.  It  faithfully  follows  the  sharp  maximum  around  the 
transition  from  regular  reflection  to  Mach  reflection  and  the  change 
of  the  character  in  the  regular  reflection  region  as  the  burst  height 
is  increased. 

Isopicnic  contours  of  peak  dynamic  pressure  are  plotted  versus 
scaled  burst  height  and  ground  range  and  compared  with  data  points  in 
Figs.  101  through  109.  The  dynamic  pressure  range  is  0.05  to  1000 
ps  i . 


SCALED  DYNAMIC  IMPULSE  VERSUS  SCALED  BURST  HEIGHT  AND  GROUND  RANGE. 

A  simple  fit  exists  for  the  integral  of  dynamic  pressure  with 
time  over  the  positive  (outward  flow)  phase,  as  defined  in  Eq.  (53). 
This  fit  is  restricted  to  the  Mach  reflection  region  and  has  the  form 


+  /  Ex  +  G 

U  'F  +  x3*61  1  +  0.22x2 


m  psi-ms  , 


(65) 


where  E  =  l83(y2  +  0.00l82)/(y2  +  0.00222), 

F  =  0.00058  exp  (9.5y)  +  0.0117  exp  (-22y), 

G  =  2.3  +  29  y/ (  1  +  1760y5)  +  25y4/(1  +  3.76y6), 

but  this  approximation  is  valid  only  for 

x  >  Xi  *  1 70i|)/ (  1  +  337^1  7 4 )  +  0.91  4ip2 , 

where  14  5  y  +  0.09, 

y  =  scaled  burst  height  H  in  kilofeet  per  cube-root  kiloton, 

=  H/IOOOW173, 

x  -  scaled  ground  range  GR  in  kilofeet  per  cube-root  kiloton, 
=  GR/IOOOW173. 
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SCALED  H08  (  KFT) 


CONTOURS-  .4.  .3.  .2.  .16  PSI 


SCALED  GR  (  KFT) 

Note:  a  =  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a,  b] ;  Smiley,  Tomayko, 
and  Ruetenik  [1982  a-d] ;  Ruetenik  [1984] ;  and  Brode  [1959b] . 


Figure  102.  HOB  isopicnic  contours  from  Kaman  AviDyne  calculations  compared 
to  fit  for  peak  dynamic  pressure  (0.15  to  0.4  psi). 


Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a,  b] ;  Smiley,  Tomayko, 
and  Ruetenik  [1982  a-d]  ;  Ruetenik  [1984]  ;  and  Brode  [1959b] . 

Figure  103.  HOB  isopicnic  contours  from  Kaman  AviDyne  calculations  compared 
to  fit  for  peak  dynamic  pressure  (0.08  to  0.15  psi). 


CONTOURS-  .38.  .07.  .06.  .06  PSI 


Note:  a  =  calculation  result. 

Source:  Smiley,  Ruetenik,  and  Tomayko  [1984a,  b] ;  Smiley,  Tomayko, 
and  Ruetenik  [1982  a-d] ;  Ruetenik  [1984] ;  and  Brode  [1959b]. 

Figure  104.  HOB  isopicnic  contours  from  Kaman  Avidyne  calculations  compared 
to  fit  for  peak  dynamic  pressure  (0.05  to  0.08  psi). 
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This  fit  differs  from  the  surface  burst  [ 1  — KT  standard  (Needham  and 
Crepeau,  1981)]  by  nearly  7  percent  at  700  to  1000  psi-ms  (y  =  0),  but 
is  within  5  percent  everywhere  else.  The  average  difference  is  less 
than  2  percent.  The  fit  is  compared  with  the  calculations  in 
Fig.  110. 

In  Fig.  110,  the  calculation  values  are  indicated  by  symbols 
(circles,  squares,  diamonds,  triangles)  and  the  fit  by  solid  curves. 
The  range  of  applicability  is  for  bursts  below  700  ft/KT1^  and  is  for 
x  >  Xt.  The  scaled  ground  range  is  then  between  300  and  5500  ft; 
i.e.,  the  scaled  impulses  are  between  2  and  2000  psi-ms/KT1 ^3 . 

DYNAMIC  PRESSURE  TIME  VERSUS  SCALED  BURST  HEIGHT  AND  RANGE. 

The  quick  fix  analytic  approximation  [Speicher,  1982]  has  not  yet 
been  replaced  by  a  more  detailed  fit.  This  fit  to  the  early  KA  cal¬ 
culations  (and  the  surface  burst)  is  cast  in  terms  of  the  overpressure 
f i t  of  Eq.  (63) • 


for  X  >  1 . 3  Xm  , 


(66) 


for  X  <  1 .3  Xm  , 


where  Q  is  the  dynamic  pressure  (horizontal  component)  in  pounds  per 
square  inch,  Xm  is  the  scaled  range  at  which  Mach  reflection  begins 
for  a  given  burst  height  (in  feet  per  cube-root  kiloton) ,  APS  is  the 
peak  overpressure  at  the  scaled  burst  height  and  range  [Eq.  (62)],  AP 
is  the  time-dependent  overpressure  at  that  position  [Eq.  (63)]  and  X 
is  the  scaled  ground  range  in  feet  per  cube-root  kiloton,  the  coeffi¬ 
cients  a  and  b  are  defined  as  2  -  2/(1  +  38l7y9)  and  2  +  1.011/(1  + 
33660y15)>  respectively,  and  y  is  the  burst  height  in  kilofeet  per 


ki loton 


One  further  correction  is  necessary  to  alter  the  effective  posi¬ 
tive  phase  duration.  That  involves  replacing  the  duration  of  the 
positive  phase  D,  from  Eq.  (63),  in  the  overpressure  fit  by 
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Figure  105.  Scaled  dynamic  impulse  fit  versus  scaled  range  and  scaled  burst  height. 


C  x  D  ms  , 


(67) 


D  = 
u 


where 


C 


89 .6y 


5.2 


4.51 


(1  +  20. 5y5,4)  (1  +  130. 7y8*6) 


Q,  2  2,1 .25 

1  2 . 8  ( x  +  y  ) _ 

,  A  2  ^  2,1.25  ’ 

1  +  3.o3(x  +  y  ) 


2.466y°‘5 

(1  +  99y2*5) 


and  x  =  X/1000  =  scaled  ground  range  in  kilofeet  per  cube-root  kiloton. 


The  peak  dynamic  pressure  and  the  dynamic  impulse  derived  from 
this  quick  fix  fit  are  not  as  accurate  as  those  given  by  Eqs.  (64)  and 
(65).  To  illustrate  the  limited  usefulness  of  this  approximation 
[Eqs.  (66)  and  (67)],  the  peak  dynamic  pressures  are  compared  with 
those  from  the  KA  calculations  in  Figs.  Ill  and  112,  all  at  the  scaled 
burst  height  of  200  ft/KT1^.  While  the  peaka  from  this  time-history 
fit  are  appreciably  low  at  the  innermost  range  plotted  in  Fig.  Ill 
(24.5  psi  versus  30  psi) ,  the  agreement  gets  better  at  larger  ranges 
(Fig.  112).  The  impulses,  shown  in  Figs.  113  and  114,  are  similarly 
poor  at  the  closest  range  shown,  but  they  are  in  good  agreement  far¬ 
ther  out  at  intermediate  ranges.  This  fit  falls  about  30  percent  low 
on  peak  and  impulse  at  the  100  psi  overpressure  range.  This  quick  fix 
time-history  fit  was  provided  as  an  analytic  expression  useful  in 
dynamic  analyses  in  a  limited  (low)  overpressure  range,  and  should  not 
be  used  at  high  overpressures.  The  expressions  for  peak  dynamic 
pressure  and  dynamic  impulse  [Eqs.  (64)  and  (65)]  are  more  accurate 
for  peaks  and  total  impulse,  but  they  do  not  provide  the  time- 
dependent  or  transient  behavior  necessary  for  calculations  of  response 
of  structures  or  vehicles.  It  is  anticipated  that  this  quick  fix 
time-history  fit  will  be  improved  in  the  near  future. 

Figure  115  shows  a  comparison  between  the  KA  calculations  and  the 
dynamic  pressure  iso-impulse  contours  from  the  integration  of 
Eq.  (66).  As  expected,  this  comparison  is  not  as  good  as  that  for 
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Dynamic  pressure  (psi) 


Dynamic  pressure  (200-ft  SHOB) 


Figure  106. 


Fit  compared  to  Kaman  AviDyne  calculation:  peak  dynamic  pressure 
for  200-ft  SHOB  versus  scaled,  close-in  ground  range. 
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Dynamic  pressure  (200-ft  SHOB) 


Figure  107-  Fit  compared  to  Kaman  AviDyne  calculation:  peak  dynamic  pressure 
for  200-ft  SHOB  versus  scaled,  intermediate  ground  range. 


1 40 


Scaled  dynamic  impulse  (psi-s/KT 


<r> 


Scaled  dynamic  impulse  (200-ft  SHOB) 


Fiqure  108.  Fit  compared  to  Kaman  AviDyne  calculation:  peak  dynamic  impulse 
for  200-ft  SHOB  versus  scaled,  close-in  ground  range. 


1 41 


Scaled  dynamic  impulse  (psi-s/KT ,/J) 


0.16 


Scaled  dynamic  impulse  (200-ft  SHOB) 


Figure  109.  Fit  compared  to  Kaman  AviDyne  calculation:  peak  dynamic  impulse 
for  200-ft  SHOB  versus  scaled,  intermediate  ground  range. 
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Figure  110.  Integration  of  Eq.  (66)  over  positive  phase  compared  to  Kaman  AviDyne 
calculations:  dynamic  impulse  versus  SHOB  and  scaled  ground  range. 


Eq.  (65),  shown  in  Fig.  110.  However,  it  corresponds  reasonably  well 
to  the  total  impulse,  while  supplying  a  time-history,  which  is 
presumably  more  useful  in  dynamic  analyses  of  structural  response. 

HEIGHT  OF  TARGET  EFFECT  ON  MAXIMUM  DYNAMIC  PRESSURE  CONTOURS. 

As  in  Figs.  101  through  109,  the  curves  of  Fig.  116  show  the 
scaled  range  to  KA  calculation  peak  dynamic  pressure  values  (at  scaled 
burst  heights  of  0,  200,  250,  400,  500,  700,  and  750  ft/KT173).  Also 
shown  in  Fig.  116  are  the  locus  of  burst  heights  and  ranges  where  the 
peak  dynamic  pressure  at  a  scaled  height  above  the  ground  of  1,  2,  3, 
5,  and  10  ft/KT1/8  is  less  than  that  at  the  ground  by  just 
10  percent.  At  points  above  the  dashed  curves  in  Fig.  116,  the  peak 
dynamic  pressures  are  presumably  even  lower  than  at  the  surface  for 
each  target  height.  Below  those  curves,  the  peak  dynamic  pressure  at 
a  target  height  above  the  surface  is  within  10  percent  of  the  surface 
values.  In  the  figure,  solid  curves  are  for  surface  values  of  peak 
dynamic  pressure.  These  curves  of  the  10  percent  limit  for  target 
heights  are  closely  fit  by  the  expression 

Rt  =  a  +  by  +  cyd/(1  +  eyf)  kft/KT1/3  ,  (68) 


where  y 
a 
b 
c 
d 
e 
f 
TH 


scaled  burst  height  in  kilofeet  per  cube-root  kiloton, 
-1.49  +  0.4073(TH)"0-81866, 

-13.1  +  1 6(TH)+0-023, 

0.06791  +  0.0363(TH)  +  3.019x1 0“8(TH)7, 

-1 .531 3  +  0.1 302 ( TH) , 

1.039(TH)  +  34.83/(TH)1 -6754, 

2.733  +  4. 1 58/ [1  +  0. 001 59(TH) 4- 246] , 
target  height  in  feet  per  cube-root  kiloton. 
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is  90  percent  of  surface  value. 


SECTION  5 

EQUATION  OF  STATE  FOR  AIR 


For  strong  blast  waves  in  the  atmosphere,  the  real  gas  (nonideal 
gas)  properties  of  air  become  important,  and  the  ideal  gas  assumptions 
frequently  prove  inadequate.  Air  molecules  enter  a  complex  energy- 
density-pressure  balance  as  they  dissociate  and  ionize  with  rising 
pressure  or  temperature.  Their  behavior  is  expressed  in  thermodynamic 
terms  from  detailed  calculations  [Gilmore,  1955,  1959;  Hilsenrath, 
Green,  and  Beckett,  1957].  The  equation  of  state  for  air  has  been 
closely  fit  by  the  author  [Brode  and  Parkin,  1963]  and  that  fit  has 
been  used  in  detailed  numerical  calculations  of  the  radiation  flow  and 
hydrodynamics  of  nuclear  bursts  [Brode,  1959a, b,  1966,  1969;  Brode 
et  al. ,  1)67],  as  well  as  in  the  KA  calculations  [Smiley,  Tomayko,  and 
Ruetenik,  1982a-d]. 

CALORIC  EQUATION  OF  STATE  FOR  AIR. 

The  usual  ideal  gas  relation  for  specific  internal  energy  can  be 
written  as 


E  = 


(69) 


in  which  E  is  the  energy  per  unit  mass,  P  is  the  pressure,  p  is  the 
density,  and  Y  is  the  ratio  of  specific  heats  (Y  =  Cp/Cv) .  Defining 
U  =  (Y  +  1)/(Y  -  1),  one  can  rewrite  Eq.  (69)  as 


2pE 

(4-1)  * 


(70) 


Using  a  dimensionless  variable  defined  as 


<t>  = 


\1  .0553 


(71  ) 
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in  which  PQ  is  the  standard  sea  level  pressure  (14.7  psi),  P  is  the 
air  pressure,  p  is  the  air  density,  p0  is  the  standard  density 
(1.293  g/1),  and  £  =  ln(p/p0)  =  lnp,  then  the  fit  is  expressible  as 

U=1  +  A  +  B  +  C  +  D  +  E  +  F  +  G  +  H+  I  +  (27<f>  +  3) /(5<f>  +  1)  ,  (72) 

where  each  literal  term  A  through  I  has  the  form 

(Mq  +  H  5  +  M  £2)<f>m(1  -  <j>  )a 

J  =  - - - n -  ‘  (73) 

(M3  +  M4^)64>  +  1 

The  coefficients  in  the  literal  terms  are  listed  in  Table  3  in  which 
powers  of  ten  are  abbreviated  as  follows:  2.236+5  means  2.236  x  105. 

This  form  is  valid  to  better  than  5  percent  almost  everywhere  for 
all  pressures  up  to  10^  psi  and  for  densities  in  the  range  10“^  < 
p/p0  <  20.  Radiation  energy  is  not  included  in  the  fit;  at  very  high 
temperatures,  the  effective  specific  heat  ratio  is  that  of  a 
monoatomic  gas  (Y  =  5/3).  At  standard  conditions,  the  ratio  is  close 
to  that  for  an  ideal  diatomic  gas  (Y  =  7/5).  Since  the  effective 
specific  heat  ratio  changes  slowly  with  pressure  or  density,  it  is 
usually  easy  to  find  the  appropriate  thermodynamic  state  by  using  a 
few  iterative  steps. 

The  data  on  which  this  fit  is  based  are  represented  by  the  curves 
of  Fig.  117;  the  corresponding  fit  is  shown  in  Fig.  118. 

THERMAL  EQUATION  OF  STATE  FOR  AIR. 

In  defining  a  temperature,  the  thermodynamic  relation  between 
pressure,  density,  and  temperature  becomes  more  complex  than  the  usual 
ideal  gas  form,  which  is: 


p  -  PSe  .  (74) 

The  gas  constant  %  becomes  a  variable.  See  Sec.  2  [Eq.  (24)  and 
Fig.  6]  for  a  fit  to  that  variability. 
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Table  3-  Caloric  equation  of  state  for  air. 
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Figure  112.  Caloric  equation  of  state  for  air--data. 
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SECTION  6 
THERMAL  RADIATION 


A  nuclear  fireball  releases  a  fair  fraction  of  the  explosion 
energy  in  a  pulse  of  light  and  heat  or  thermal  radiation.  This  burst 
of  radiant  energy  is  complex  and  variable  in  time,  spectrum,  and 
geometry.  It  is  coupled  with  the  nuclear  radiation  and  the  blast 
wave.  It  has  at  least  two  major  peaks,  and  varies  considerably  with 
yield,  altitude,  and  atmospheric  conditions.  When  close  to  the 
ground,  the  earth  material  ingested  into  the  fireball  markedly  affects 
the  thermal  radiation.  When  burst  at  high  altitude,  the  thiner  atmos¬ 
phere  leads  to  more  rapid  radiation  as  inner  hot  regions  are  exposed 
earlier,  yet  the  greater  transparency  leads  to  low  emissivities  and 
low  thermal  power  (radiation)  at  higher  temperatures  or  earlier  in  the 
fireball  expansion  and  cooling. 


AIRBURST 


For  nuclear  explosions  in  the  atmosphere,  but  not  on  the  ground, 
the  thermal  radiation  may  be  characterized  by  a  double  pulse,  the 
first  of  which  is  so  short  that  it  contains  less  than  half  a  percent 
of  the  total  energy  released,  yet  its  peak  power  is  comparable  to  that 
in  the  second  pulse  which  contains  more  than  a  third  of  the  total 
yield.  The  major  or  second  pulse  is  well  approximated  by  a  simple 
time  dependence: 

Thermal  Power 


P 


2m  a 


(75) 
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where  P*  =  power  or  rate  of  fireball  heat  release, 
rr  =  tr/t2max> 

tr  =  time  after  burst,  in  seconds, 

t^o  =  time  to  second  maximum  in  power  radiated  (airburst) 


0  44  t 

«  0.0417‘W  s,  ±11% T 


(76) 


P  =  power  at  second  maximum  from  airbursts, 
ma 

=  3.18  W0-56  KT/s,  ±  34%  (77) 

W  =  yield,  in  kilotons. 


Equations  (76)  and  (77)  are  from  Glasstone  and  Dolan  [1977],  and 
represent  an  empirical  approximation  to  observed  and  calculated  times 
to  maximum  and  power  at  maximum. 

The  pulse  approximated  by  Eq .  (75)  and  the  partial  integral  of 
that  pulse  over  time  are  illustrated  in  Fig.  119.  The  integral  of  the 
power  fit  to  a  time  tpQma  =  10  tma  includes  91  percent  of  the  integral 
to  very  late  times,  i.e.,  it  leaves  only  9  percent  to  be  radiated 
after  that  time.  The  curve  with  which  it  is  compared  [Glasstone  and 
Dolan,  1977]  leaves  20  percent  beyond  t^Qma-  The  curve  plotted  for 
the  fit  (Fig.  119)  has  been  renormalized  to  80  percent  at  tqOma- 
The  expression  for  thermal  power  [ Eq .  (75)],  when  integrated 
analytically,  leads  to: 


E  ( r ) 


—  arctan  (/2r/(l  -  r^) 

J2 


t ^  4-  fir  +  1 

r2  -  J2r  +  1 


(78) 


t The  ±11%  indie  ates  that  a  band  around  the  empirical  values  of  ±11 
span  90%  of  the  atmospheric  test  data. 
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Figure  114.  Fit  compared  to  data  average:  normalized  thermal  radiation  power 
pulse  and  thermal  energy  radiated  versus  normalized  time. 


This  function  (normalized  at  tioma  to  Percent  of  the  total 
thermal  emission)  is  plotted  as  the  fit  in  Fig.  119. 

First  Peak 

The  main  burst  of  thermal  energy  is  generally  preceded  by  a  sharp 
pulse  of  light,  as  the  fireball  shock  breaks  through  the  air  surround¬ 
ing  the  bomb  that  has  been  altered  (ionized)  by  the  early  nuclear  and 
X-ray  radiation  and  begins  to  cool  as  it  expands  further.  That  first 
pulse  contains  less  than  1/2  percent  of  the  total  yield,  and  peaks  in 
a  time  measured  in  milliseconds.  The  time  to  the  first  peak  is 
proportional  to  the  cube-root  of  the  yield  (compatible  with 
hydrodynamic  scaling): 


t 


fma  — 


ms . 


(79) 


Thermal  Minimum 

A  minimum  in  radiated  power  occurs  between  the  first  and  main 
peaks  when  the  expanding  and  cooling  shock  front  becomes  sufficiently 
transparent  that  the  hotter  interior  air  of  the  fireball  is  exposed, 
allowing  an  increase  in  the  rate  of  radiant  energy  emission.  This 
time  to  minimum  power  for  air  bursts  depends  on  optical  properties 
(mean  free  paths  for  light  passing  through  air)  as  well  as  fireball 
dimensions.  An  approximate  empirical  scaling  from  Kieth  and  Sachs 
[1985]  is  as  follows: 


t  . 
min 


3.13-W 


o.44 


ms. 


±12% 


(80) 


Thermal  Energy  Fraction 

For  bursts  low  in  the  atmosphere,  the  radiated  energy  at  a  time 
ten  times  the  time  to  maximum  can  be  approximated  by 
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£  * 

f  10ma  '  E  1 0ma  /W  =  °-277  KT  ±  19  %  (81) 

If  radiation  at  times  beyond  ten  times  time  to  maximum  is  included  in 
the  thermal  fraction,  the  expression  becomes  [Glasstone  and  Dolan, 
1977] 


#  £ 

f  «,  -  E  „  /W  =  0.35.  (82) 


An  earlier  fit  to  the  thermal  energy  (in  visible  light)  as  a 
function  of  yield  resulted  in  the  expression  [Brode,  1968]: 

f*^  —  0.33  +  0.0038-W1  /2/( 1  +  0.031 6-W1 /2)  .  (83) 


This  latter  form  (Eq.  83)  predicts  a  rise  in  the  fraction  with  in¬ 
creasing  yield  (from  33  to  42  percent  between  1  and  10,000  KT) ,  while 
Glasstone  and  Dolan  [1977]  hold  the  fraction  fixed  at  35  percent.  In 
contrast,  the  recent  revision  of  the  Thermal  Radiation  section 
(Chapter  6)  of  the  Defense  Nuclear  Agency  manual  compiled  by  Kieth  and 
Sachs  [1985]  predicts  a  decrease  with  increasing  yield  (although  they 
state  that  f*  =  0.33  is  equally  appropriate.  An  analytic  fit  to  their 
curves  for  a  low  airburst  is: 

#  ?  ii 47 

f  -  34. 9/(  1  +  0 . 001 62  •  X  *  )  (84) 

00 


in  which  X  =  ln(W)  (W  in  KT). 

These  three  predictions  are  illustrated  in  Fig.  120,  showing  a 
divergence  between  the  values  as  yield  is  increased.  Although  a  good 
deal  of  scatter  exists  in  the  atmospheric  test  data,  they  do  not 
appear  to  support  the  reduction  of  thermal  fraction  with  increasing 
yield  predicted  by  Kieth  and  Sachs  [1985]. 
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Figure  115.  Thermal  fraction  versus  yield,  near  sea  level  airburst  predictions. 


Radiant  Exposure 


The  intensity  or  fluence  of  fireball  thermal  radiation  received 
at  a  distance  from  an  airburst  falls  off  as  the  square  of  the  dis¬ 
tance,  and  is  also  attenuated  by  absorption  and  scattering  in  the 
intervening  air.  It  may  be  affected  by  reflections  from  ground  cover 
(e.g.,  snow)  or  clouds  or  by  absorption  in  intervening  cloud  layers, 
and  will  depend  on  burst  height  and  yield.  From  an  airburst  weapon, 

the  thermal  fluence  at  a  distance  can  be  expressed  in  the  following 
f  orm: 


* 

Q 


B 


T  /[  4u 


*  2 

(r  T] 


cal/cm 


(85) 


in  which 


attenuation  through  clouds 

enhancement  due  to  snow  or  cloud  reflections 
total  radiated  energy  (cal) 
transmissivity  of  the  atmosphere 
distance  from  burst  (cm) 


£ 

The  factor  A  ,  representing  the  attenuation  through  clouds  can 
range  from  a  value  of  unity  (for  no  clouds)  to  as  much  as  1/10  for 
heavy  cloud  cover  between  burst  and  target.  [1/10  <A*  <1] 

The  enhancement  due  to  snow  cover  can  be  as  great  as  1.9,  and 
nearly  as  great  for  reflections  from  a  cloud  deck  above  the  burst. 
[1  <E*  <1.9] 


It  is  conventional  to  express  the  energy  release  in  kilotons; 
there  are  lO1^  calories  in  a  kiloton.  The  fraction  of  the  energy 
released  that  appears  in  thermal  radiation  for  a  low  altitude  airburst 
[f*  E«/W]  is  given  variously  by  Eqs.  82  through  84. 

The  atmospheric  transmission  is  a  complex  function  of  the  radia¬ 
tive  properties  of  the  air  molecules  and  the  suspended  contaminants, 
but  can  be  usefully  approximated  without  reference  to  its  spectral 
dependence  by  an  absorptive  (extinction)  property  and  a  forward  scat¬ 
tering  contribution.  The  absorption  can  be  characterized  by  a  term 
[exp  (-a*  •  R*/V*)],  and  the  scattering  factor  as  (1  +  B*  •  R*/V*) ,  in 
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which  V*  is  a  visibility  or  atmospheric  visual  distance.  Thus, 

T*  =  ( 1  +  B**R*/V*) *exp(-a**R*/V*)  /  (86) 

. 

with  R*  and  V*  in  the  same  length  units  (e.g.,  statute  miles). 

With  ground  range  in  statute  miles,  Eq .  (85)  becomes 

*  *  #  x  p  p 

q  =  3 • 07 • A  -B  -f  *T  *W/(R  )  cal/cm  ,  (87) 

where  R^j_  =  [(H*)^  +  (GR*)2]1/2  mi,  | 

H^j_  =  Burst  height  mi, 

GR*  =  Ground  Range  mi. 

With  distances  in  kilometers,  the  expression  becomes 

£  $###  #p  p 

Q  =  7. 96 • A  -B  -f  -T  -W/(R  ,  )  cal/cm  .  (88) 

km 


For  nautical  miles  (used  in  some  targeting  work): 

Q*  =  2.32*A**B**f**T**W/(R*  )2  cal/cm2  .  (89) 

nm 

Ranges  are  given  in  kilofeet  elsewhere  in  this  compendium  of  fits;  the 
thermal  fluence  versus  slant  range  in  kilofeet  is  written: 

Q*  =  85.7*A**B*‘f**T**W/ (R*  )2  cal/cm2  (90) 


4-  %  £ 

The  constants  a  and  8  depend  on  the  atmospheric  moisture  and 
particulate  contents,  but  best  values  are  ct*o^2.9,  8*^1. 9.  A 
reasonable  range  for  each  might  be  1.8<a<3.2  and  1.25<B<2.1. 
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SURFACE  BURST 

The  thermal  Dulse  from  a  surface  burst  is  less  regular  and 
generally  less  intense  than  that  from  an  airburst.  Soil  and  other 
surface  material  drawn  into  the  fireball,  intervening  topography  and 
structures,  and  greater  distortion  of  the  fireball  radiating  surface 
due  to  reflections  and  precursors  often  degrade  the  thermal  pulse  from 
surface  and  near-surface  bursts.  The  surface  burst  power  pulse  may 
have  the  same  general  shape  as  the  air  burst  [Eq.  (75)],  but  the  time 
to  maximum,  the  power  at  maximum,  and  the  total  energy  emitted  are 
different.  [Kieth  and  Sachs,  1985]: 


*  ,  0.464 

t  ~  0.041 9  *W  s  , 

ma 


*  0  56 

P  ~  1  .35-W  KT/s  , 

ma 


*  *  n  0 . 02  4 

f  =  E  .  _  /W  ~  0 . 1  1  8  -W 

1 0ma  —  1  Oma 


*  *  ,0.024 

f  =  E  /W  ~  0.  1 49-W  KT 

o°s  °°s 


(91) 

(92) 

(93) 

(94) 


Graphical  data  from  Fig.  6II-5  from  Kieth  and  Sachs  [1985]  is  fit 
by  the  following  expression  when  extrapolated  to  zero  burst  height: 

f*  ~  0.04475  +  0.0021  61  *  A 3/ ( 1 +0.01575*X3) ,  (95) 

COS 


in  which,  again,  X  =  ln(W). 

In  addition,  the  time  to  first  maximum  is  less  regular,  and  the 
time  to  minimum  is  greater: 


0.44 

t  ~  4 . 25  *W  tns  • 

min 


(96) 
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Thermal  Partition — Transition  between  Surface  and  Airburst 

Kieth  and  Sachs  [1985]  give  some  guidance  as  to  the  transition 
between  surface  and  airburst  thermal  partition  in  their  Table  6II-2. 
When  fit  with  an  analytical  expression,  this  transition  may  be 
presented  as: 


f*  ~  f*  +  (f*  -f*  )*(h/4)2/[1  +  ( h/4) 2 ]  . 

t  s  s  s 


(97) 


In  this  formula,  h  is  the  scaled  burst  height  in  meters  per  cube  root 
kiloton,  fs  is  the  surface  burst  partition  [Eq.  (95)],  and  f  is  the 
airburst  thermal  partition  [Eq.  (83)].  This  transition  is  illustrated 
as  a  function  of  scaled  burst  height  for  yields  from  1  KT  to  10  MT  in 
Fig.  121 . 

BURIED  BURSTS 

Burial  quickly  quenches  the  thermal  radiation  from  a  nuclear 
burst,  but  at  shallow  depths  of  burst,  a  fraction  of  the  energy  still 
radiates  from  the  contaminated  fireball.  Based  on  the  limited  ex- 
pereince  of  a  f ew  near  surface  and  shallow  buried  bursts,  the  follow¬ 
ing  fit  approximates  the  thermal  energy  partition  for  such  shallow- 
buried  bursts: 


f*  ~  . 35 • { 1 -1 /[I +( (h+5.925)/8.1 44)^]} , 


for  -4.89  ^h  <100m/(KT) 


1/3 


f  ~  0  for  h  ^-4.89  . 
u 


(98) 


HIGH  ALTITUDE  BURSTS 


For  bursts  above  about  14,00-ft  altitude 
of  decreased  air  density  on  fireball  behavior 
atmospheric  density  (p)  is  often  expressed  in 


(=4.3  km),  the  influence 
becomes  appreciable.  The 
terms  of  the  standard 
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Figure  116.  Approximate  transition  from  surface  burst  to  airburst  thermal  radiation 
fraction  versus  burst  height  for  yields  from  1  KT  to  10  MT  Eq.f(97)l. 


sea-level  atmospheric  density  ( p~1 .225  grams  per  liter),  and  can  be 
approximated  [U.S.  Standard,  1962]  by  the  function: 

n  =  p/p  ~  [1-*-0.5405-(h/1  0)°*9852/(1  +0. 15^5-(h/1  O)3'393) 
o 

3  909  17  74 

+  5.591  *(h/100r  /( 1  +36- 32  -  (h/1  00)  )] 

•exp(-0 . 1 4856 *h)  ,  (99) 


in  which  h  is  the  altitude  in  kilometers.  This  expression  fits  the 
tabular  data  to  better  than  3  percent  up  to  30  km.  It  is  good  to 
better  than  6  percent  up  to  115  km. 

The  following  expressions  approximate  the  thermal  pulse  charac¬ 
teristics  in  the  atmosphere  (h<30  km): 

*  0  4S  0  2? 

t  .  =  4.3-W  -n  J  ms  ,  (100) 

min 

*  0  45  0  90 

t  =  40-W  -ri  ms  ,  (101) 

ma 

P*  =3.1 -8*55  /n*54  KTs  .  (102) 

ma 

The  time  to  first  maximum  is  nearly  independent  of  altitude  of 
burst,  so  that  Eq.  (79)  can  be  used  for  any  height.  The  fraction 
emitted  as  thermal  radiation  can  be  roughly  approximated  by  [Keith  and 
Sachs,  1985]: 


* 

1  Oma 


~  0.276/n 


0.034 

f 


(103) 


f 


# 

00 


2*  0.35/n 


0.34 


(104) 
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This  ratio  of  thermal  energy  radiated  to  the  total  yield  of  a 
nuclear  explosion  at  altitude  has  been  approximated  earlier  [Brode, 
1968],  based  on  high  altitude  tests  and  a  f ew  numerical  calculations, 
for  densities  as  low  as  n  ~  0.001  and  for  yields  from  1  KT  to  1 0  MT 
as: 


f  ^0.27  +  0.06«n  +  0.0038-W1 /2/(  1 +0.031  6‘W1 /2) 

+  1 1 4 • n/ ( 1 +82000-n2)  .  (105) 

The  expression  in  Eq.  (105)  indicates  an  increase  in  thermal 
radiation  fraction  with  increasing  yield,  as  the  emissivity  of  the 
fireball  remains  high  longer  and  to  lower  fireball  temperatures.  That 
allows  more  radiation  to  be  emitted.  The  fraction  shows  a  slight 
initial  drop  with  atmospheric  density,  because  at  lower  densities,  the 
fireball  forms  a  shock  later  (thereby  becoming  a  strong  radiation 
source  later),  and  yet  it  becomes  transparent  earlier.  Eventually, 
radiative  expansion  in  a  rarefied  atmosphere  creates  such  a  large 
fireball  that  it  can  again  radiate  effectively  in  the  visible.  At 
extreme  altitudes,  the  lack  of  atmosphere  does  not  allow  for  reradia¬ 
tion  in  the  visible,  and  the  thermal  fraction  once  more  decreases. 
Figure  122  shows  the  fit  [Eq.  (105)]  as  a  function  of  yield  for  burst 
altitudes  up  to  30  kilometers. 

The  thermal  fraction  versus  yield  and  altitude  of  burst,  as 
depicted  by  Glasstone  and  Dolan  [1977]  can  be  approximated  by  the 
form : 
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-  35  +  a.h9‘5/(l+b-h10)  +  c-h3  , 


a  =  9.621  *1  0~9-[X4/(  1  +2.0737*1  0  6-A9) 


+  3. 1 88 - X 2/ ( 1  +0.009269-A2) ]  , 


+  4 . 857 • X / ( 1  +0.2895-A)  ], 


c  =  0.0002042-(  1  -0.0776-A  +  0 . 0001 1 81  . x 3 *  882 


(106) 


where  A  =  ln(W),  and  h  =  altitude  in  kilometers. 

This  approximation  [Eq.  (106)]  is  plotted  in  Fig.  123  versus 
yield  for  burst  heights  up  to  30  km. 

The  thermal  fraction  versus  yield  and  altitude  as  depicted  by 
Kieth  and  Sachs  [1985]  is  at  variance  with  these  previous  estimates 
[Eqs.  (105)  and  (106)].  It  can  be  fit  by  the  form: 


f 


2  3 

a  +  b  *h  +  c  *h  +  d  *h 


(107) 


00 


a  =  34.86/(1 +0.001618-A 


2.  447 


)  , 


in  which 


b  =  0.1125/(1 +0.08972-A 


1.545 


-  0.005319  , 


c  =  0.01861 -0.009506’A+0.005402-A2-0.001008-A3 
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Figure  1 1 8 .  Thermal  radiation  fraction  of  total  yield  versus  yield  for  burst  altitudes 
ud  to  30  KM  [E1JW,  19771. 


again  with  \  =  ln(W)  (W  in  KT)  and  h  =  altitude  of  burst  in 
kilometers.  This  fit  is  illustrated  in  Fig.  124. 

The  three  predictions  approximated  in  Eqs .  (105),  (106)  and  (107) 
are  based  on  many  of  the  same  test  data  and  numerical  calculations,  so 
the  differences  are  not  easily  justified.  The  test  data,  as  well  as 
most  of  the  calculations,  contain  considerable  scatter  and 
variability,  and  the  information  at  altitude  is  sparse.  However,  more 
and  better  calculations  are  possible,  and  could  help  to  further 
resolve  the  uncertainties.  The  three  preditions  are  plotted  versus 
altitude  for  1  KT,  100  KT,  and  10  MT  in  Figs.  125,  126,  and  127, 
respectively.  All  the  forms  lead  to  an  increase  in  the  thermal  frac¬ 
tion  with  altitude  at  high  altitudes  (h>l4km),  but  the  Kieth  and  Sachs 
[1985]  values  decrease  with  yield.  The  test  data,  although  limited, 
do  not  seem  to  support  a  decreasing  thermal  fraction  with  increasing 
yield  at  any  altitude. 
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Figure  120.  Comparison  of  thermal  fractions  as  Dredicted  by  three  formulae 
(Eas.  (101),  (102),  (103))  versus  burst  altitude  for  1  KT. 
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Figure  121  .  Comparison  of  thermal  fractions  as  predicted  by  three  formulae 
(Eqs .  (101),  (102),  (103))  versus  burst  altitude  for  100  KT. 
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Figure  122  .  Comparison  of  thermal  fractions  as  predicted  by  three  formulae 
(Eqs.  (101),  (102),  (103))  versus  burst  altitude  for  10  MT. 
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APPENDIX  A 

OVERPRESSURE  VALUES  AS  FUNCTION  OF 
GROUND  RANGE,  BURST  HEIGHT,  AND  TIME 


X 

Y 

a  -  t 

AP 

18.2 

0.0 

0. 042 

1 07926.6 

28.1 

0.0 

0.363 

10113.7 

37.9 

0.0 

1  .22 

2069.0 

47.1 

0.0 

0.0045 

27696.0 

66.3 

0.0 

6.39 

306.5 

80.3 

0.0 

19.8 

78.34 

0.0 

25.0 

0.137 

55580.0 

19.7 

25.0 

25.3 

9.284 

32.8 

25.0 

0.00417 

255273.1 

49.2 

25.0 

0.046 

1 0462.5 

64.7 

25.0 

0.189 

1 5346. 0 

83.0 

25.0 

57.8 

7.799 

97.6 

25.0 

124.0 

1  .265 

0.0 

50.0 

0.041 6 

T1 0106.0 

6.72 

50.0 

0.116 

43739.3 

24.2 

50.0 

3.07 

270.9 

41  .5 

50.0 

56.1 

0.1174 

55.9 

50.0 

0.137 

21 674.8 

72.6 

50.0 

0.251 

1 51 68.9 

90.6 

50.0 

1 1  .7 

70.46 

14.9 

22.5 

0.01 39 

201704.4 

33.1 

18.9 

1  .35 

1 397.0 

37.7 

8.22 

0.21  4 

1 6365.0 

82.8 

1 41 0.0 

6.95 

8.374 

157.0 

257.0 

0.448 

375.8 

393.0 

682.0 

10.9 

24.23 

1 250.0 

220.0 

22.1 

6.593 

1 1 60.0 

2520.0 

0.538 

3.1 185 

31 40.0 

4260.0 

13.4 

1  .233 

NOTE: 


X  =  scaled  ground  range  (ft/KT 
Y  =  scaled  burst  height  (ft/KT17^), 
a  -  x  -  scaled  time  -  time  of  arrival  ( 
AP  =  overpressure  (psi). 


ms/KT17^) 


f 
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APPENDIX  B 

FORTRAN  PROGRAM  FOR  OVERPRESSURE 
VERSUS  TIME,  BURST  HEIGHT,  AND  GROUND  RANGE 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  PT< Y ,X , SI GMA , DELTAP) 


THIS  IS  A  FORTRAN  IMPLEMENTATION  OF  THE  ANALYTIC  EXPRESSION  FOR 
PRESSURE  TIME  HISTORY  < BRODE  AND  SPEICHER,  MAY  1986) 


THE  PARAMETERS  ARE: 


X 

Y 

SIMGA 


RANGE 


SCALED  GROUND 
SCALED  HOB 
SCALED  TIME 
**  NOTE  THAT  SIGMA 
**  ANALYTIC 
**  SIGMA  IN 
**  WHICH  IS 


< FT/KT  **< 1/3) ) 
<FT/KT**< 1/3) ) 
(MSEC/KT  **< 1/3) ) 
IS  THE  TIME  AFTER 


BURST  IN  THE 
EXPRESSION.  FOR  PURPOSES  OF  CALCULATION 
THIS  ROUTINE  IS  TIME  AFTER  TIME  OF  ARRIVAL, 
THEN  ADDED  TO  THE  CALCULATED  TIME  OF  ARRIVAL 


** 

** 

*•* 


DELTAP 


PRESSURE 


<PSI  ) 


IF  PRESSURE  IS  DESIRED  IN  KPA  THEN  MULTIPLY  RESULT 
BY  100/14.504 

NOTE  THAT  LIMITS  ARE  PLACED  ON  VALUES  SUCH  AS  X,Y,Z,  ETC.. 
THIS  IS  DONE  TO  AVOID  OVERFLOWS  AND  THE  VALUES  ARE  MACHINE 
DEPENDENT . 


IF  <X.LT. IE-9)  X=*l . E-9 
IF  < Y . LT . 1 E-9)  Y»1 .E-9 
CAPR  =»  <X*X  +  Y*Y)**.5 


R  =  CAPR/ 1000 

Z  =  Y/X 

IF  (Z.GT.100.)  Z5*  100 

CALL  PPEAK<X , Y , CAPR , Z , DELTPS) 

XM  =  1 70 *Y  /  <1.  +  60 .*<Y**.25) )  +  2 . 89*< < Y/l 00 ) **2 . 5) 

U  =*  <0.543  -  21 . 8*R  +  386*<R**2)  +  2383*< R**3) ) *< R**8) 

#  /  <  2 . 99E-1 4  -  1 .91E-10*<R**2)  +  1 . 032E-6*< R**4)  - 

II  4.43E-6*<R**6)  +  <1.028  +  2.087»R  +  2.69*<R**2))  *  <R**8)) 
TAU  =  U 

IF  < X . LT . XM )  GOTO  100 

W  =  <1.086  -  34 . 605*R  +  486.3*<R**2)  +  2383*<R**3)) 

#  *  <R**8)  /  <3.01 37E-1 3  -  1 . 21 28E-9*< R**2) 

♦♦  +  4 . 1  28E-6*<  R**4)  -  1  . 1  1 6E-5*<  R**6)  +  <1.632  + 

#  2 . 629*R  +  2 . 69*<  R**2) )  *  <R**8)) 

TAU  =*  U*XM/X  +  W*<  1  ~  XM/X) 

100  CONTINUE 

SIGMA  »  SIGMA  +  TAU 

52  =  1  -  1 5 . 1 8*< < Y/l 00 ) **3. 5)  /  <1  +  1 5 . 1 8*< < Y/l 00 ) **3 . 5) )  - 


179 


#  < 0 . 02441 *< <Y/1 . E6) **2)  /  <1  +  ?000*< <Y/1 00 ) ) **7) )  *  Cl . E 1 0  / 

#  (0.441  +  < <X/100)**10) ) ) 

CAPD  =  <<1<640700  +  24629*TAU  +  41 6 . 1 5*<TAU**2)  )  /  (10880  + 

#  619.76*TAU  +  <  TAU**2) ) )  *  <0.4  +  0 . 00 1 204*<TAU**1  . 5)  / 

#  <1  +  0 . 001559*<TAU**1 .5) )  +  <0.6126  +  0 . 5486*<TAU»* . 25)  / 

#  <1  +  0.00 357 * < TAU** 1.5))  -  3.47*<TAU**0 . 637)  /  <1  + 

#  5.696*<TAU**0 .645) ) )  *  S2> 

S  =  1  -  1 100*< <Y/100>**7)  /  <1  +  1 100*< <Y/100)**7) )  - 

#  <  2 . 441 E-l 4*Y*Y  /  <1  +  9000*< < Y/l 00 ) **7) ) )  *  (1.E10  /  <0.441  + 

#  < <X/100)**10) ) ) 

F2  =  <0.445  -  5 . 44*( R**l .02)  /  <1  +  1 00000*< R**5 . 84) )  + 

#  7 . 571 *<2**7. 15)  /  <1  -  5. 135*<Z**12.9) )  -  8 . 07*< Z**7 . 31 )  / 

#  <  1  ♦  5.583*<Z**12.23) ) )  *  < 0 . 435*< < Y/l 0 ) **1  . 26)  /  <1  + 

#  0 .03096*< <Y/10)**3. 12) ) )  *  <1  -  0 .00001?*<TAU**8)  /  <1  + 

#  0 .00001 9*<TAU**8) ) ) 

F  =  <0 .01477*<TAU**.75)  /  <1  +  0 . 005836*TAU)  +  7 . 402E-5*< TAU**2 . 5) 
tt  /  <1  +  1 .429E-8*<TAU**4.75) )  -  0.216)  *  S  +  0.7076  - 
tt  3.077E-5*<TAU**3)  /  <1  +  4 . 367E-5*<TAU**3) )  +  F2  -  <0.452  - 

#  9.94E-7*<X**4. 13)  /  <1  +  2 . 1 868E-6*<X**4 .13)))  *  <1  - 

#  1 .5397E-4*<Y**4.3)  /  <1  +  1 . 5397E-4*< Y**4 . 3) ) ) 

G  =  10  +  <77.58  -  64.99*<TAU»*. 125)  /  <1  +  0 . 04348*<TAU**5) ) )  *  S 
H  =  3.003  +  0 . 0560 1 *TAU  /  <1  +  1 . 473E-9*<TAU**5) )  +  <0.01769*TAU  / 

#  <1  +  3 . 207E-1 0*<TAU**4 . 25) )  -  0 . 0320?*<TAU**1 . 25)  /  <1  + 

#  9.91 4E-8*<TAU**4) )  -  1.6)  *  S  -  0 . 1 966*<TAU**1 . 22)  /  <1  + 

#  0 .767*<TAU**1 .22) ) 

B  a  <F  *  < <TAU/SIGMA)**G)  +  <1  -  F)  *  < (TAU/SI GMA) **H) )  * 
tt  <1  -  (SIGMA  -  TAU)/CAPD) 

IF  (X.LT.XM  .OR.  Y.GT.380.)  GOTO  1000 
XE  =*  3 . 039*Y  /  <  1  +  0.0067*Y) 

AK  =»  ABS<  <X  -  XM)/(XE  -  XM)  ) 

IF  (AK.GT.50.)  AK=50 

D2  =»  2.99  +  31  240* <  < Y/l  00  )  **9 . 86)  /  <  1  +  1  5530*<  ( Y/l  00  )  **9  .  87)  ) 

D  =*  0.23  +  0 . 583*Y*Y  /  <26667  +  Y*Y)  +  0.27*AK  +  <0.5  -  0.583*Y*Y  / 

#  <26667  +  Y*Y)  )  *  <AK**02) 

A  =  <D  -  1)  *  <1  -  < AK**20 )  /  <1  +  ( AK**20 ) ) ) 

AJ  »  11860  *  (SIGMA  -  TAU)  /  <Y  *  < <X  -  XM)**1.25)) 

IF  (AJ.GT.200)  A J=20 0 
U  =  1  +  ( 0 .0037 44* <  (Y/10)**5. 185)  / 

#  <1  +  0 .004684*< <Y/10)**4. 189) ) 

#  +  0 .004755*< (Y/10)**8.049)  /  <1  +  0 . 003444*< < Y/l 0 ) **7 . 497) ) 

#  -  0 . 04852*<  < Y/l 0 ) **3 .423)  /  <  1  +  0 . 03038*< < Y/l 0 ) **2 . 538) ) )  * 

4  < AJ**3)  /  <6.13  +  <AJ**3) )  *  <1  /  (1  +  9 . 23*< AK**2) ) ) 

C3  =  1  +  < 1 .094*<AK**.738)  /  <1  +  3 . 687*< AK**2 . 63) )  *  <1  - 

#  83 . 0 1 *<  < Y/l 00 ) **6 . 5)  /  <1  +  1 72 . 3*<  < Y/l 00 ) **6 . 04) )  -  0.15))  * 

#  <1  /  <1  +  0 .5089*(AK**13) ) ) 

C2  »  2300 0*<  < Y/l 00 ) **9)  /  <  1  +  23000 *<  < Y/l 00 ) **9) ) 

TEMP  =  <XX 100) **4 
C  -  <1.04  -  0 . 02409*TEMP  /  <1  + 

4  0 . 0 23 1 7 *TEMP ) )  *  <AJ**7)  /  <<1  +  A)  *  <1  + 

#  0 . 923*<AJ**8. 5) ) )  *  <C2  +  <1  -  C2)  *  <1  - 

#  0 . 09*<  AK**2 . 5)  /  <  1  +  0 . 09*(  AK**2 . 5)  )  )  ) 

#  *  C3  *  < 1  -  << (SIGMA  -  TAU) /CAPO) **8) ) 

DELTAP  =  DELTPS  *  <1  +  A)  *  <  B*L>  +  C) 

RETURN 

1000  CONTINUE 

DELTAP  =  DELTPS  *  B 

RETURN 

END 
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SUBROUTINE  PPEAK(X ,Y , CAPR , Z , DELTPS) 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CAPR  =»  (X*X  +  Y*Y)**(l/2) 
Z  =  Y/X 

DELTPS  =  PEAK  PRESSURE 


THIS  IS  A  FORTRAN  IMPLEMENTATION  OF  THE  ANALYTIC  EXPRESSION  FOR 
PEAK  OVERPRESSURE  ( BRODE  AND  SPEICHER,  MAY  1986) 


THE  PARAMETERS  ARE: 


X  =*  SCALED  GR 
Y  =  SCALED  HOB 


(PSI ) 


<  FT/KT **( 1/3) ) 

<  FT/KT **( 1/3) ) 
(  FT/KT **(  1/3)  ) 


R  =  CAPR/ 1 000 

A  =»  1.22  -  (3.908*Z*Z)  /  (1  +  810. 2*Z**5) 

B  =  2.321  +  <  6 . 1  95*(  2**1 8)/( 1  +  1 . 1 1 3*(Z**1 8) ) )  - 

#  (0 .03831*<Z**17)/( 1  +  0 .02415*<Z**17) ) )  +  (0.6692/(1  + 

#  41 64*( Z**8) ) ) 

BB  =  . 0629*( (X/Y) **8 . 34)  /  (1  +  . 00509*( (X/Y) **1 3 . 05) )  * 

#  . 05*Y  /  (1  +  2 . 56E-8*( Y **5) ) 

C  =  4.153  -  ( 1 . 1 49*(Z**1 8) )/  (1  +  1 . 641 *( Z**l 8) )  -  (1.1  /  (1  + 

#  2 . 771 *( Z**2 . 5) ) ) 

D  =  (-4.166)  +  25.76*(Z**1 .75)  /  (1  +  1 . 382*( Z**l 8) )  +  3.257*2  / 

#  ( 1  +  3.21 9*Z ) 

E  =  1  -  0 . 004642*(Z**1 8)  /  +  0  P>«9986*(Z**1 3)  ) 

F  =  0.6096  +  2 . 879* ( Z**9 . 25)/  (1  +  2 . 359*( Z**l 4 . 5) )  -  17.15*Z*Z  / 
«  ( 1  +71 . 66*( Z**3) ) 

G  =»  1.33  +  5 . 361  *Z*Z  /  (1  +  .  31  39*(  Z**6) ) 

H  =  -(0.2905  +  64 . 67*( Z**5) )  /  (1  +  441 . 5*( Z**5) )  -  1.389*2  /(I 

#  +  49 . 03*( 2**5) )  +  8.908*(Z**1 .5)  /  (1  +  1 54 . 5*( 2**3 . 5) )  + 

M  1 .094*(CAPR**2)  /  ((0.7813E9  -  1.234E5*CAPR  +  1 20 1 *( CAPR**1 . 5)  + 

#  ( CAPR**2) )  *  (1  +  2*Y) ) 

P  =  1 . 8008E-7*( Y**4)  /  (1  +  0 . 0002863*(Y**4) )  -  2.121*Y*Y  / 
(794300+  ( Y **4 . 3) ) 

Q  =  5.18  +  8 . 864*( Y**3 . 5)  /  ( 3 . 788E6  +  (Y**4) ) 

OELTPS  31  (  1 0  •  47)/(  R**A)  +  (B  -  BB)/(  R**C)  +  (D*E)/(1  + 

F*(R**G))  +  H  +  P/(R**Q) 


RETURN 

END 
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APPENDIX  C 
LIST  OF  SYMBOLS 


air  density,  in  kilograms  per  cubic  meter; 

Eqs.  (6)  through  (8) 

air  temperature,  in  degrees  kelvin;  Eq.  (24) 

ambient  air  density 

ambient  air  sound  speed 
=  1 .0872  kft/s 

ambient  air  specific  heat  ratio 
*  1.400 

ambient  air  temperature  =  273.2  K; 

Eqs.  (20),  (21) 

ambient  pressure,  in  pounds  per  square  inch 
»  14.7  psi  at  sea  level 

burst  height,  in  feet 

duration  of  negative  overpressure  phase, 
in  milliseconds 

duration  of  positive  overpressure  phase, 
in  milliseconds;  Eqs.  (43)  through  (48) 

duration  of  positive  dynamic  pressure,  or 
outward  flow;  Eq.  (52) 

dynamic  pressure,  in  pounds  per  square  inch 
(pu2/2);  Eqs.  (15)  through  ( 1 8 ) ,  (56),  (57), 
(63),  (64),  (66) 

fireball  temperature  maximum,  in  103°C;  Eq.  (59) 

free-air  or  surface-burst  range,  in  kilofeet; 
Eqs.  (36),  (40),  (41),  (55) 

ground  range,  or  shock  radius,  in  feet 


P 

9 

po 

2o 

^o 


H 


d; 


K 


Q 


sr 

GR 


NOTE:  Subscripts  "o"  refer  to  ambient  air  (preshock)  conditions; 
subscripts  "s"  refer  to  shock  conditions;  and  subscripts  "R"  refer 
to  reflected  shock  values. 
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impulse  in  dynamic  pressure  positive  phase, 


in  pounds  per  square  inch  millisecond; 

Eqs .  (53)  through  (55),  (65) 

ij 

impulse  in  positive  overpressure  phase, 
in  pounds  per  square  inch  millisecond; 

Eqs.  (47),  (48) 

rP 

interface  between  regular  and  Mach 

reflectioo*  approximate;  Eq.  (64)  Xq 


locus  of  scaled  burst  heights  and  scaled  ranges 
for  second  peak  equal  to  shock  overpressure, 
in  kilofeet  per  cube-root  kiloton;  Eq.  (63) 

Xe 

lower  limit  to  range  of  validity  of  dynamic 
impulse  fit;  Eq.  (65) 

*i 

Mach  number  (u/C) 

M 

Mach  reflection  onset,  in  kilofeet 
per  cube-root  kiloton;  Eq.  (63) 

Xm 

overpressure,  in  millions  of  pounds  per 
square  inch 

C 

overpressure,  in  pounds  per  square  inch; 

Eqs.  (49),  (50),  (63) 

AP 

overpressure,  in  thousands  of  pounds  per  square 
square  inch  =  APS/1000 

IT 

particle  velocity 

U 

peak  density 

ps 

peak  dynamic  pressure  at  transition  between 
regular  and  Mach  reflection;  Eq.  (64) 

Qm 

peak  dynamic  pressure,  in  pounds  per  square  inch; 

Eqs.  (16)  through  (18) 

Qs 

peak  overpressure,  in  pounds  per  square  inch; 

Eqs.  (33)  through  (35),  (37),  (38),  (50),  (62) 

APS 

peak  particle  velocity,  in  kilofeet  per  second; 

Eqs.  (12)  through  (14) 

us 

NOTE:  Subscripts  "o"  refer  to  ambient  air  (preshock) 

conditions ; 

subscripts  "s"  refer  to  shock  conditions;  and  subscripts  "R"  refer 
to  reflected  shock  values. 
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peak  pressure,  in  pounds  per  square  inch; 

Eqs.  (2)  through  (6) 

positive  dynamic  pressure  duration;  Eq.  (58) 

positive  overpressure  duration;  Eqs.  (50),  (51) 

power,  or  rate,  of  heat  release 

pressure  =  AP  +  P0,  in  pounds  per  square 
inch 

radiant  fluence,  in  calories  per  square 
centimeter 

range  for  height  of  target  at  which  dynamic 
pressure  is  lower  by  10  percent  from  surface 
value  for  corresponding  burst  height;  Eq.  (67) 

range  for  thermal  radiation,  in  statute  miles 

ratio  of  scaled  burst  height  to  scaled  ground 
range  (y/x);  Eqs.  (62),  (63) 

reflected  peak  overpressure  (normal  reflection), 
in  pounds  per  square  inch;  Eq.  (30) 

reflected  pressure,  in  pounds  per  square 
inch;  Eqs.  (26)  through  (30) 

reflected  shock  velocity;  Eq.  (25),  (26) 

reflection  factor  =  APR/AP3;  Eqs.  (30) 
through  (32) 

scaled  burst  height,  in  feet  per  cube-root 
kiloton  =  H/m 

scaled  burst  height,  in  kilofeet  per  cube-root 
kiloton  =  Y/1000 

scaled  ground  range,  in  feet  per  cube-root 
kiloton  =  GR/m 

scaled  ground  range,  in  kilofeet  per  cube-root 
kiloton  =  X/1000 


P 


s 


D' 


D 

P* 

P 


Rrp 

Rr 


z 


APr 

PR 

Ur 

RF 

Y 

y 

x 


x 


NOTE:  Subscripts  "o"  refer  to  ambient  air  (preshock)  conditions, 

subscripts  "s"  refer  to  shock  conditions;  and  subscripts  "R"  refer 
to  reflected  shock  values. 
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scaled  slant  range,  in  feet  per  cube-root 
kiloton  =  (X2  +  Y2)1/2 

R 

scaled  slant  range,  in  kilofeet  per  cube-root 
kiloton  =  (x2  +  y2)1/2;  Eq.  (33) 

r 

scaled  time  after  detonation  (t/m),  in  milli¬ 
seconds  per  cube-root  kiloton 

a 

scaled  time  of  arrival  (T/m),  in  milli¬ 
seconds  per  cube-root  kiloton 

T 

shock  density  ratio 

^S 

shocked  air  specific  heat  ratio 

^s 

shock  temperature 

6S 

shock  temperature  increase,  in  degrees 
centigrade;  Eqs.  (20)  through  (23) 

A0S 

shock  velocity,  in  kilofeet  per  second; 

Eqs.  (10),  (11) 

us 

specific  heat  at  constant  pressure 

cp 

specific  heat  at  constant  volume 

cv 

specific  heat  ratio  =  Cp/Cv;  Eq.  (9) 

Y 

specific  internal  energy,  in  ergs  per  gram; 

Eqs.  (4),  (5),  (69),  (70) 

E 

specific  volume,  in  cubic  centimeters  per  gram 

-  1/p 

V 

speed  of  sound 

c 

standard  deviation 

E 

thermal  gas  constant  =  PV/0 


thermal  gas  constant  for  air  at  shock  conditions 

s 

thermal  gas  constant  for  air  at  ambient  conditions 

0 

NOTE:  Subscripts  "o"  refer  to  ambient  air  (preshock) 

conditions ; 

subscripts  "s"  refer  to  shock  conditions;  and  subscripts  "R"  refer 
to  reflected  shock  values. 
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V 


thermal  gas  constant  ratio  =  /  =  pv/(  e)- 

Eqs.  (24),  (25)  0  ’ 

thermodynamic  variable,  equation  of  state 
for  air,  dimensionless;  Eq.  (71) 

time  after  detonation,  in  milliseconds 

time  after  detonation,  in  seconds  (thermal) 

time  of  arrival,  in  milliseconds;  Eqs.  (39) 
through  (41 ) 

time  to  maximum  fireball  temperature;  Eq.  (60) 

time  to  negative  phase 

total  energy  radiated 

transmissivity 

visibility,  in  statute  miles 

yield,  in  kilotons 

yield  scaling  factor  =  cube-root  yield  in 
cube-root  kilotons 


<P 

t 

t  r 


W 


m 


a  = 

-0.086  - 
9n  *1  o 

’4;  Eq.  (25) 

6  = 

6(o, n) ;  Eq. 

(25) 

n  = 

O 

CL 

\ 

Q. 

u  - 

i 

>- 

+ 

>■ 

1  )  =  1  +  2Ep/P;  Eq 

to  = 

(t  -  t)/d*= 

(a  -  t)/D;  Eq.  (56) 

=  h  +  0.09  (a  parameter  in  Qg  and  1^  fits  versus  HOB); 
Eqs.  (64),  (65) 

S  =  Inn 


NOTE:  Subscripts  "o"  refer  to  ambient  air  (preshock)  conditions; 

subscripts  "s"  refer  to  shock  conditions;  and  subscripts  "R"  refer 
to  reflected  shock  values. 
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